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1.  Introduction 


Analysis  of  conditions  for  crack  initiation  and  fracture  toughness  in  the 
case  of  metallic  materials  under  various  loading  rates  are  the  problems  of 
great  interest  in  fracture  mechanics.  For  Mode  I  and  Mode  III  several 

investigations  have  been  carried  out  to  reveal  the  loading  rate  effect  (  Ki, 

km)  on  fracture  toughness  (ATjc,  ATnic)-  A.t  the  same  time,  systematic 
analyses  have  been  performed  for  majority  of  structural  materials  to  study 

strain  rate  (e)  effect  on  the  yield  stress  {Gy),  ultimate  stress  (Ou)»  material 
strain  hardening,  etc.  Using  of  such  data,  some  correlations  between 
fracture  toughness  and  parameters  of  stress-strain  diagram  at  different 
strain  rates  have  been  derived,  (Klepaczko,  1990).  The  theoretical 

background  of  such  correlations  consists  of  solutions  for  elastic  and  elastic 
plastic  stress-strain  fields  near  a  crack  tip,  such  as  Williams'  (1957)  linear- 
elastic  solution  for  crack  tip  fields  at  a  stationary  crack;  approximate 
analyses  of  Irwin  (1960),  Dugdale  (1960)  (often  refered  to  as  the 
Barenblatt-Dugdale-Bilby-Cottrell-Swinden  model),  Duffy  (1969)  for  the 
size  and  shape  of  the  crack  tip  plastic  zone;  solutions  of  Tuba  (1966)  as 
well  as  Hutchinson  (1968a)  and  Rice  and  Rosengren  (1968)  (refered  to  as 
the  HRR-solution)  for  crack  tip  fields  in  elastic-plastic  materials.  Analyses 
of  stress  distribution  at  the  tip  of  a  stationary  crack  under  dynamic  loading 
have  been  performed  by  Freund  (1990)  and  Lee  and  Freund  (1990). 

Application  of  the  HRR-solution  and  criterion  of  cleavage  fracture  in 
the  form  (Ritchie  et  al.,  1973) 

iCj  =  when  =  Op  at  .r  =  /p 

results  in  the  following  relationship 


^ic  “ 


J<+l)/2 

■  a  I(N)  /p  ’ 

(i-v2)j:^-i 

|l/2 


(11) 


where  Gp  is  the  critical  cleavage  stress;  /p  is  characteristic  distance  ahead  of 
the  crack  tip;  cr^^is  the  stress  component  in  Canesian  coordinates  {x,y) 
with  the  origin  at  the  crack  tip  (Fig.l);  a  and  N  are  material  constants  in 
the  Ramberg-Osgood  equation 


—  =  —  ot  — 


(1.2) 


with  E=g^/E-,  E  and  v  are  Young's  modulus  and  Poisson's  ratio;  Z  and  I(N) 

Y  Y  =  yy 
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are  coordinate  function  and  numerical  factor  in  the  HRR-solution  (see 
Chapter  3). 

Another  important  application  of  fracture  mechanics  arises  from  a 
relationship  between  loading  rate  and  specimen  thickness  for  fracture 

toughness  measurements,  i.e.  B=f(^i),  The  reason  for  limitations  on  the  B 
value  turns  out  from  requirements  of  small  scale  yielding  or,  equivalently, 
restrictions  on  the  crack  tip  plastic  zone  size,  r^.  Thus,  making  use  of 

Irwin's  (1960)  correction  for  planj  strain,  one  finds 


^  SJT  (Oy  j 


(1.3) 


The  well-known  ASTM  limiting  condition 

B>  2.5  (1.4) 

is  found  to  be  equivalent  to  the  following  one 

B>24rp  (1.5) 

Since  the  tp  value  is  diminishing  with  increasing  the  loading  rate, 
specimens  smaller  in  dimensions  can  be  employed  in  tests  at  high 

loading  rates.  An  example  of  correlation  between  B  and/Cj  is  given  by 
Klepaczko  (1990). 

It  should  be  noted  that  the  majority  of  experimental  and  theoretical 
results  on  rate  sensitivity  of  fracture  toughness  are  obtained  for  cracks  in 
Mode  I  and  Mode  III  (see  Klepaczko,  1990).  For  in-plane  shear  c.racks 
(Mode  II)  such  solutions  are  not  available.  No  much  progress  is  achieved  in 
experimental  techniques  for  measurements  of  fracture  toughness  (A^iic»  .^iic) 
in  Mode  II.  Only  a  few  works  are  known  on  experimental  studies  of  Mode 
II  crack  behavior  at  quasi-static  case  (Banks-Sills  and  Arcan,  1986;  Banks- 
Sills  and  Sherman,  1991;  Davies  et  ai,  1985;  Hoyniak  and  Conway,  1979; 
Maccagno  and  Knott,  1992;  Mall  and  Mol,  1991;  Tohgo  et  ai,  1989;  Tohgo 
and  Ishii,  1992)  and  dynamic  case  (Kalthoff  and  Winkler,  1987;  Mason  et 
al.,  1992).  Some  important  features  of  crack  behavior  in  Mode  II  which 
follow  from  those  works  are  outlined  below: 

(1)  Difficulties  in  creating  pure  Mode  II  conditions.  Only  in  quasi¬ 
static  loading,  pure  Mode  II  can  be  expected  to  occur  for  some  specimens 
and  load  schemes.  More  often,  mixed  Mode  I-II  crack  tip  deformation 
(^^{>0)  or  crack  surfaces  contact  (A(i<0)  are  observed.  In  the  latter  case, 
application  of  the  /-integral  and  HRR-solution  for  characterization  of  crack 
tip  fields  becomes  questionable. 

(2)  Different  character  of  fracture  initiation  for  specimens  of  different 
materials  and  with  different  load  schemes.  Thus,  some  authors  observe 
crack  initiation  in  direction  of  about  70. ..80°  with  respect  to  the  initial 
crack  line,  which  is  usually  refered  to  as  brittle  fracture  mechanism  (Sih, 


5 


1973;  1974);  in  this  case  the  crack  initiation  angle  can  be  predicted  by 
criteria  of  maximum  tensile  stress  or  minimum  strain  energy  density.  In 
the  majority  of  investigations  (Banks-Sills  and  Sherman,  1991;  Maccagno 
and  Knott,  1992;  Tohgo  and  Ishii,  1992)  crack  propagates  in  direction  close 

to  the  original  crack  (or  notch)  plane,  which  is  refered  to  as  ductile  fracture 

mechanism. 

(3)  Transition  in  fracture  mechanism  with  increasing  loading  rates.  It 
has  been  observed  by  Kalthoff  and  Winkler  (1987)  in  dynamic  tests  of 
steel  that  after  certain  limit  of  the  loading  rate  the  change  in  the  crack 
initiation  angle  (from  70°  with  respect  to  the  ligament  at  lower  loading  rate 
to  about  10°  at  higher  loading  rate)  occured.  This  transition  has  been 
supposed  to  occure  due  to  formation  of  localized  shear  bands. 

Also  some  difficulties  arise  in  numerical  simulation  of  crack  behavior 
at  high  loading  rates.  An  analytic  solution  for  the  stress  field  at  a 

stationary  crack  tip  in  Mode  I  and  II  under  dynamic  loading  has  been 

derived  by  Lee  and  Freund  (1990)  for  the  case  of  linear-elastic  material 
and  idealized  geometry  of  a  cracked  body,  i.e.  an  edge  crack  in  a  half- 
plane,  Similar  analysis  for  specimens  of  finite  size  and  elastic-plastic 
material  seams  to  be  very  complicated  with  prime  attention  paid  to  use 
adequate  stress-strain,  strain-rate,  temperature  constitutive  relations 
governing  material  behavior  at  the  crack  tip  in  high  strain  rates  (see 
Klepaczko,  1987). 

As  an  alternative  to  exact  solution  of  elastic-plastic  stress-strain 
fields  at  a  crack  tip  under  dynamic  loading,  an  approximate  quasi-static 
approach  can  be  applied  as  follows.  Equations  (1.1), (1.3), (1.4),  as  well  as 
other  solutions  for  evaluation  of  Kjc,  tp,  B  and  stress  and  strain 
components,  involve  material  parameters  (aY,a,N)  which  are  in  general 
rate  sensitive.  Since  those  material  parameters  can  be  found  for  certain 
material  as  functions  of  strain  rate,  their  introduction  into  equations  for 
fracture  mechanics  parameters  leads  to  approximate  assessments  of 
material  behavior  at  a  crack  tip  for  different  loading  rates.  It  is  obvious 
that  such  a  quasi-static  analysis  takes  no  account  of  inertia  effect, 
temperature  and  strain  rate  redistribution  due  to  presence  of  a  crack.  In 
fact,  when  an  exact  analysis  cannot  be  avoided,  more  general  constitutive 
relations  are  to  be  used  along  with  numerical  methods. 

In  the  present  report,  the  quasi-static  analysis  is  carried  out  to  study 
strain  rate  effect  on  elastic-plastic  stress-strain  fields  at  a  crack  tip  in 
Mode  II  for  materials  with  different  rate  sensitivity;  mild  steel,  titanium 
alloy  and  aluminum  alloy.  Evolution  of  the  crack  tip  plastic  zone  size  and 
shape  with  increasing  the  strain  rate  from  10’^^  to  10^  s'^  is  analysed.  Also 
a  brief  review  of  solutions  on  elastic  and  elastic-plastic  crack  tip  fields  in 
Mode  II  and  results  of  Kjic  or  Jhq  measurements  are  presented. 
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2.  Characterization  of  Stress  and  Strain  Fields  and  Plastic 
Zone  at  the  Crack  Tip  in  Mode  II  Via  Small  Scale 
Yielding  Approximation 

The  small  scale  yielding  (SSY)  concept  is  succesfully  applied  to  fracture 
mechanics  for  approximate  analysis  of  plasticity  effects  on  crack  behavior. 
If  the  SSY  concept  is  considered,  a  small  plastic  zone  of  radius  r=rp  is 
supposed  to  occur  at  the  crack  tip  with  a  stress  field  within  and  near  the 
plastic  zone  (r  <  rp+5)  determined  from  the  general  elastic-plastic  solution, 
while  at  the  distance  r  >  rp+5  from  the  crack  tip  the  stress  components  are 
governed  by  the  asymptotic  linear-elastic  representation  (Williams,  1957; 
Irwin,  1958) 

(0)  +  0(r®)  (2.1) 

^!2■xr 

Here  are  the  stress  components  in  Cartesian  coordinates  (x.y  )  with  the 

origin  at  the  crack  tip  (Fig.l);  (r,0)  are  the  polar  coordinates  with  the  origin 
at  the  crack  tip  (Fig.l);  Xj^^^are  specified  functions  of  the  polar  angle;  index 
"p"  is  refered  to  the  crack  tip  deformation  mode,  p  =  I,II,III  (I  -  opening 
mode,  II  -  in-plane  shear  mode.  III  -  anti-plane  deformation);  is  the 
stress  intensity  factor  for  corresponding  deformation  mode.  The  second 
term  in  the  right  side  of  eq.(2.1)  is  limited  at  r->0. 

Evidently,  the  necessary  condition  for  validity  of  the  SSY 
approximation  is  the  requirement  of  small  size  of  the  plastic  zone  in 
comparison  to  the  crack  length  a,  body  thickness  B,  ligament  size  H=W-a 
(W  being  the  width  of  a  cracked  section),  i.e. 

tp  « <3  ,  tp  «  B  ,  tp  «  H 

Practically,  in  measuring  the  plane  strain  fracture  toughness,  the  conditions 
discussed  above  are  satisfied  with  the  ASTM  requirement 

B  >  2.5 

which  corresponds  to  rpi0-04B,  when  Irwin's  plasticity  correction  is  taken. 

The  use  of  the  SSY  approximation  allows  for  linear-elastic  fracture 
mechanics  methods  and  criteria  (e.g.,  Arp=#rpc)  be  extended  to  analyse  the 
elastic-plastic  material  behavior.  In  particular,  within  the  SSY 
approximation  some  efficient  solutions  for  the  crack  tip  plastic  zone  size 
and  shape  have  been  derived  (Irwin,  1960;  Dugdale,  1960;  Duffy  et  al., 
1969;  Larsson  and  Carlsson,  1973;  Rice,  1974).  Some  approximate 
assessments  of  the  crack  tip  plastic  zone  in  Mode  II  are  discussed  below 
with  prime  attention  paid  to  the  change  of  tp  due  to  variation  of  the 
applied  stress  level,  tappi/Ty,  which  makes  it  possible  to  incorporate  the 
rate  sensitivity  of  the  yield  stress,  Ty.  into  analysis  which  follows. 


2.1.  Mode  II  Crack  Tip  Stress  Field  in  Small  Scale  Yielding 


Asymptotic  terms  of  the  Imear-elastic  stress  field  at  the  crack  tip  in  Mode 
II  are  given  by  (Williams,  1957;  Irwin.  1957;  Broek.  1987;  Kanninen  and 
Popelar,  1985) 

Oij  =  x„(e) 

ylKT 

In  Cartesian  coordinates  one  has 

V  ■  6  r-,  9  30  , 

Lxx  =  -  sm—  (2  +  cos—  cos —  ) , 

2  2  2 

v  -0  0  30 

Zyy  =  Sin-  cos-  cos—  ,  (2.2) 

=  cos—  { 1  -  sin—  sin — )  , 

2  2  2 

Z;;  =  v(Z^^  +  Zyy)  for  plane  strain. 

=  0  for  plane  stress, 

T  _  y  _  A  . 

~  “it  “  » 


in  cylindrical  coordinates 

Sr:  =  sin—  (3ccs’-  -  2)  , 

2  2 

lee  =  -  3sin|-  cos^^  ,  (2.3) 

Sr0  =  cos^  (3cos2^  -  2)  , 

7  7 

Zzz  =  v(Srr  +  See)  for  plane  strain, 

Szz  =  0  for  plane  stress, 

Se^j  =  S^r  =  0 

the  principal  stresses  are 


<^1,2=  (-sin— 

±  7/ 1  -(3/4)sin'0  )  , 

V  2jir  “ 

03  =  -2v  sin^ 

for  plane  strain. 

03=0 

for  pleine  stress 

The  stress  intensity  factor,  ATjj  ,  is  a  measure  of  the  amplitude  of  the  crack 
tip  stress  field.  In  the  presence  of  plasticity,  the  7- integral  (Cherepanov, 
1979;  Rice,  1968a)  is  found  to  be  an  adequate  parameter  characterizing  the 
crack  tip  field  amplitude  (Hutchinson,  1968a).  According  to  Rice  {1968a), 
the  7- integral  is  represented  by 
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7=  [(Wdy  -  CTijnj^)  d5  (2.5) 

j 

r 

where  f  is  an  arbitrary  contour  surrounding  the  crack  tip  with  the  origin 
and  end  at  lower  and  upper  crack  surfaces  respectively,  nj  are  components 
of  the  outer  normal  to  F,  u;  are  displacement  components  {x\=^x,  X2=yy,  W  is 
the  strain  energy  density, 


W  =  jaijdeij 
0 


(2.6) 


Under  conditions  of  traction-free  crack  surfaces  and  absence  of  body 
forces,  the  7-integral  is  path-independent.  This  property  is  employed  to 
derive  several  fracture  mechanics  solutions  both  in  the  elastic  and  elastic- 
plastic  cases.  The  7-integral  was  introduced  into  the  elastic-plastic  crack  tip 
field  solutions  (Hutchinson,  1968a;  Rice  and  Rosengren,  1968)  and,  thus, 
used  as  a  fracture  criterion  when  plasticity  induced  effects  are  significant. 

If  the  SSY  approximation  is  assumed,  the  7-integral  and  stress 
intensity  factor  for  Mode  II  are  related  by 


•^ii 


l_v2  2 


for  plane  strain. 


for  plane  stress. 


(2.7) 


2.2.  Irwin's  Plastic  Zone  Correction 


Consider  crack  of  length  a  subjected  to  in-plane  shear  by  stresses  '^appi 
infinity  (Fig. 2).  The  presence  of  the  plastic  zone  of  size  rp  ahead  of  the 
crack  tip  can  be  taken  into  account,  according  to  Irwin  (1960),  by 
increasing  compliance  of  the  cracked  body.  To  introduce  plasticity,  an 
effective  crack  of  length  aeff=‘i  +  5i  ts  considered  (Fig. 2)  with  the  stress 
intensity  factor  equal  to  Following  the  SSY  approximation,  the  stress 

at  a  distance  .r>62  ahead  of  the  effective  crack  tip  is  expressed  by 

T  _  ^Il.eff 

Ijy  -  - 

y2itx 

The  distance  62  is  found  by  equating  x^y  =  Fy  (with  Xy  be  the  yield  stress  in 
shear) 

5  =±(^ILe£f)2 
Ztc  Xy  ' 

Equilibrium  condition  (equality  of  areas  5,  and  5,  in  Fig.2)  leads  to 


=  2-\  /  ^  ^»  efr  _  i 

\  271  Xy  7t  Xy  ^ 
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The  SSY  requirement 
finally  find 


rp  «  a  is  equivalent  to  thus,  one  can 


= 


i 

:c 


(2.8) 


Formula  (2.8)  may  be  applied  to  approximate  assessment  of  the 

plastic  zone  size  ahead  of  the  crack  tip  in  Mode  II  both  in  plane  stress  and 

plane  strain.  Indeed,  analysis  of  plastic  constraint  factor  (p.c.f.)  (Broek. 
1987)  results  in 

p.c.f.  =  i./\  3 

for  both  plane  stress  and  plane  strain  conditions,  expecting  rp(6=0)  be 
independent  of  the  stress  state.  Note  that  this  feature  of  the  crack  tip 
plastic  zone  behavior  in  Mode  II  differes  substantially  from  that  in  .Mode  I. 
For  the  latter  an  effective  yield  stress  in  plane  strain  is  three  times  of  that 

in  plane  stress  (Broek,  1987).  Further  analysis  of  Irwin  (1960)  has 

revealed  the  ratio  of  p.c.f.  values  for  plane  strain  and  plane  stress  be  equal 


4^ 

to  \8  which  leads  to  the  following  approximate  assessments 

in  plane  stress. 

K  Oy 


Tg  =  in  plane  strain. 

3  It  ay 

.A.S  a  result,  the  area  of  plastic  zone  for  Mode  I  in  plane  strain  is  one  order 
smaller  that  in  plane  stress.  In  contrast,  plane  stress/plane  strain 
transition  does  not  affect  significantly  the  size  and  area  of  the  crack  tip 
plastic  zone  in  .Mode  II. 


2.3.  The  Dugdale  Model 

This  model  has  been  originally  applied  to  approximate  evaluation  of  the 
plastic  zone  size  ahead  of  a  Mode  I  crack  tip  in  plane  stress  (Dugdale, 
I960).  With  notes  made  above,  th-'re  is  no  significant  difference  in  plastic 
zone  size  for  Mode  II  crack  in  plane  stress  and  plane  strain. 

.A.S  in  Irwin's  approach,  the  effective  crack  length  aejf=ti+rp  is 
considered  under  actual  load  or  under  stresses  ^appi  applied  to  crack 
surfaces,  as  it  follows  from  the  superposition  principle  (Bueckner,  1958). 
At  the  surfaces  of  the  effective  crack,  in  small  vicinity  of  its  tip,  0  <  r  <  rp  , 
compressive  stresses  -CTy  (^ot-  opening  mode  crack)  or  shear  stresses  -ly 
(for  in-plane  shear  mode)  are  applied  to  restrict  respectively  crack  opening 
or  sliding.  The  length  tp  which  is  associated  with  the  plastic  zone  size  is 
determined  from  the  condition  of  no  singularity,  i.e.  K^=Q.  For  a  crack  in 
Mode  II  it  implies  that 


^Ob  ~  0  , 


1  0 


where  K\i^  and  are  stress  intensity  factors  for  the  effective  crack 
under  actual  stresses  t^p-i  (index  ”a")  and  stresses  -Xy  acting  within  the 
distance  tp  at  the  ctack  tip  (index  ”b"). 

Making  use  of  ’'cference  solutions  for  stress  intensity  factor  (Tada  et 
aL.  1973;  Savruk.  ’988),  one  can  find  rp for  some  crack  geometries  in  Mode 
11: 


(1)  Through  crack  of  length  2a  in  an  infinite  plate: 

=  sec  f-  -  1 


12  xv 


) 


(2)  Semi-infinite  crack  in  an  infinite  plate: 

rp  _  (Xappl^Xy 

a 


(2.9) 


(2.10) 


1  '(tappl^^Y 

Here  the  plastic  zone  size  is  normalized  by  the  crack  length,  which  is  finite 
for  a  physical  crack. 

(3)  Edge  crack  of  length  a  in  a  semi- infinite  plate. 

In  this  case  the  size  of  plastic  zone  depends  on  the  applied  load  by  the 
following  relationship 

= _ =_ _  d  ^ 


XY 


I.12I5X  a+tr 


(2.11) 


with 


1 


I(^)  =  J(l-Ti2)-i'2{i+(i_Ti2){0.2945-0.39l2Ti2+0.7685Ti^-0.9942TiM).5094TiS)}dTi 


In  deriving  eq.(2.ll)  Green’s  function  for  ATn  obtained  by  Hartranft  and  Sih 
(1973)  is  used. 

Note  an  important  difference  in  results  following  from  Irwin’s  and 
Dugdale's  models.  With  the  dimensionless  parameter  Rp  defined  as 


the  Irwin  solution  yields 


K 


with  constant  Rp-value  for  any  load  level,  Xappi/xy,  and  independent  of 
crack  geometry.  In  Fig. 3  comparison  is  made  of  Dugdale's  and  Irwin's 
solutions  for  Rp  at  a  tip  of  a  semi-infinite  crack  in  an  infmite  plane  and  an 
edge  crack  in  a  half-plane.  One  can  see  that  Dugdale's  solution  strongly 
depends  on  the  stress  level,  while  Rp-values  for  different  crack  geometries 
nearly  coincide.  More  realistic  prediction  of  the  plastic  zone  size  based  on 
Dugdale's  model,  comparing  to  Irwin's  solution,  has  been  noticed  by  Broek 
(1987)  for  the  Mode  1  conditions. 

Both  Irwin's  and  Dugdale's  approaches  allow  for  determination  of  tp 
only  ahead  of  the  crack  tip.  To  find  the  shape  of  plastic  zone,  i.c.  tp  as  a 
function  of  the  polar  angle  0,  the  linear-elastic  stress  field  in  the  form  of 


1 1 


solution  (2.2)-(2.4)  and  appropriate  yield  criterion  can  be  employed  within 
the  SSY  approximation.  The  yield  criterion  is  usually  expressed  as  the 
Huber-Mises  condition. 

(01-02)- +  (02-03)- -r  (G3-C1)- =  2a^  with  Oy=\  3  ty  .  (2.12) 

or  the  Tresca  condition, 

-‘aia.x  =  {  |cii-cr2l  .  IO2-O3I  ,  l03-0!j  }  =  Oy  with  Oy=2ty  (2.13) 

In  what  follows,  some  approximate  solutions  for  the  plastic  zone  at  the 
crack  tip  in  Mode  II  are  considered  which  are  based  on  yield  criteria  (2.12) 
and  (2.13). 


2.4.  Plastic  Zone  Geometry  in  Mode  II  Based  on 
the  Asymptotic  Linear- Elastic  Stress  Field 


Consider  a  crack  of  Mode  II  in  plane  strain.  .Making  use  of  the  Huber-Mises 
yield  condition  (2.12)  and  expressions  (2.4)  for  principal  stresses,  one  can 
find  the  following  equation  governing  the  plastic  zone  geometry  near  the 
crack  tip 


rp(9)  = 


_L  [  3 

6n.  ly  j  L 


-  z9 
+  s  1  n  — 

T 


O  ^  ^  1 

-  sin-^  -  4v(l-v)sin“  — 

4  2  ! 


(2.14) 


If  the  y-integral  instead  of  the  stress  intensity  factor  is  employed  to 
characterize  the  crack  tip  stress  field,  then  eq.(2.14)  takes  the  form 


rn(0)  = 


1 


JuE 


6a(l-v')  x3 


j  +  sin  — 


9  .  ^ 

—  Sin  d 
4 


4v(  1  -v)sin 


(2.14a) 


with  yji  and  /fn  being  related  by  eq.(2.7). 

When  the  Tresca  yield  criterion  (2.13)  is  used,  it  is  necessary  to  note 
that  according  to  eqs  (2.4) 


0;<0,  .  02<03; 

at  the  same  time  one  of  the  conditions 

Oj  >  03  or  Oj  <  03 

must  hold  depending  upon  the  values  of  v  and  0.  Thus,  the  yielding  occures 
if  one  of  the  following  conditions  is  satisfied: 

i  (Oj  -  Oo)  =  Xy  or  ^  (03  -  O2)  =  Xy 

These  result  in 


with 


rp  =  max  Upi;  rp2} 


(2.15) 
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fpi  - 


In  , 


J(1  sin^e) 


fpi  - 


Sti 


'■II 


Vl^^(3/4)sin^0~~ +  (l-2v)sin^ 


or,  in  terms  of  the  7- integral, 

r  -  ^ 

Tpi  - 


*^II^  <•  1  3  „;_2  , 


,  (l-^sm0) 

271(1 -V*)  x3  4 


(2.15a) 


fpi  - 


1  -^ii^ 

8jc(1-v')  Xy 


V 1  -(3/4)sin’e  +(l-2v)sin^ 


The  results  of  calculations  of  the  plastic  zone  shape  in  plane  strain 
based  on  the  Huber-Mises  and  the  Tresca  yield  criteria,  eqs  (2.14)  and 
(2.15),  are  presented  in  Fig. 4  for  Poisson's  ratio  ranging  from  0  to  C.5. 
(Since  the  plastic  zone  is  symmetrical  with  respect  to  the  crack  line,  only 
its  variation  at  0  <  6  <  tc  is  analysed  throughout  this  report.)  The  Huber- 
Mises  yield  criterion  predicts  (Fig. 4a)  a  smaller  plastic  zone  size  ahead  of 
the  crack  than  behind  the  tip,  with  the  condition  rp(e=O)=rp(0=7t)  valid  only 
at  vs:0.5.  One  can  conclude  from  Fig.4a  that  only  slight  differences  in  the 
plastic  zone  shape  ahead  of  the  crack  tip  exist  with  variation  of  Poisson's 
ratio;  also,  if  v>0.3  the  plastic  zone  is  nearly  symmetrical  with  respect  to 
the  y-axis  (9=jt/2).  The  Tresca  yield  criterion  predicts  (Fig. 4b)  symmetrical 
plastic  zone  shape  with  respect  to  the  y-axis  for  all  v -values  greater  than 
0.3;  this  case  corresponds  in  Fig. 4b  to  the  curve  with  v=0.33.  If  v<0.3,  only 
slight  dependence  of  the  plastic  zone  shape  on  Poisson's  ratio  in  the  range 
of  polar  angle  from  about  27c/5  to  n  is  observed. 

In  accordance  to  eqs  (2.4)  for  principal  stresses,  the  plastic  zone 
geometry  in  plane  stress  can  be  obtained  from  eqs  (2.14)-(2.15a)  putting 
Poisson's  ratio  equal  v=0.  Thus,  the  curves  in  Fig.4a,b  corresponding  to  v=0 
can  be  considered  at  the  same  time  as  plane  stress  solutions. 

Since  the  stress  intensity  factor,  ATjj,  is  proportional  to  the  applied 
stress,  Xappi*  plastic  zone  size  in  formulae  (2.14)  and  (2.15)  is 
proportional  to  (Xappi/xv)^.  Consequently,  decreasing  the  stress  level  (due  to 
increase  of  Xy  with  Xappi=const  or  decrease  of  x^ppj  with  XY=const)  leads  to 

diminishing  the  plastic  zone  proportionally  to  (Xappi/Xy)  ,  without  any 
changes  in  its  shape.  This  conclusion  is  in  contradiction  to  the  finite- 
element  results  by  Banks-Sills  and  Sherman  (1990)  which  demonstrate 
changes  in  the  shape  of  plastic  zone  with  increasing  load  level;  namely,  a 
significant  decrease  in  the  ratio  of  rp(0=7c)  to  rp(0=O)  with  increase  of  the 
load  level  has  been  observed. 
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2.5.  Influence  of  Non-Singular  Terms  in  the  Elastic  Crack  Tip 
Stress  Field  on  the  Plastic  Zone  Geometry 

The  SSY  approximation,  as  discussed  at  the  beginning  of  this  chapter,  is 
valid  only  in  a  limited  range  of  the  level  of  applied  load,  tappi/xY-  Larsson 
and  Carlsson  (1973)  and  Rice  (1974)  have  found  that  neglecting  non- 
singular  terms  in  the  series  expansion  (2.1)  of  the  stress  field  near  a  crack 
tip  leads  to  significant  errors  in  evaluation  of  the  plastic  zone  size.  For 
Mode  I,  a  two-parameter  crack  tip  field  characterization  has  been 
suggested  in  those  papers 

Oif  (r.e)  =  r  tD(e)  +  Ty  w)  . 

NIkt 

with  =  T  =*=  0,  Try  =  Tyy  =  0  (Williams,  1957;  Irwin,  1957;  Kanninen  and 
Popelar,  1985).  The  T-effect  has  appeared  to  be  different  for  various 
geometries  of  cracked  specimens  resulting  in  different  configurations  of 
the  plastic  zone  (Larsson  and  Carlsson,  1973).  At  the  same  time,  no  T-effect 
on  the  y-integral  has  been  found  by  Rice  (1974).  Another  important  piece 
of  information  which  follows  from  Rice's  (1974)  analysis  is  that,  when  T- 
effect  being  taken  into  account,  deviation  of  the  plastic  zone  size  from  the 
SSY  solution  of  Irwin  (1960)  is  proportional  to  Oappi/cty  .  while  the 
difference  is  proportional  to  (o^pp/oy)  ^  with  T-stresses  being  neglected 
(e.g.,  Dugdale's  model).  Thus,  the  two-parameter  solution  of  Larsson  and 
Carlsson  (1973)  and  Rice  (1974)  predicts  more  rapid  increase  of  the  plastic 
zone  size  with  the  load  level  than  previous  models.  With  T-effect  taken 
into  account,  the  boundary  layer  formulation  has  been  modified  and  some 
important  corrections  to  elastic  and  elastic-plastic  crack  tip  fields  have 
been  obtained  for  both  pure  Mode  I  and  mixed  Mode  I-II  conditions 
(Larsson  and  Carlsson, 1973;  Du  and  Hancock,  1991;  Betegon  and  Hancock, 
1991;  Du  et  al.,  1991). 

To  study  the  effect  of  non-singular  stress  components  on  plastic  zone 
behavior  in  Mode  II,  let  us  represent  the  crack  tip  stress  field  as  follows 
(Williams,  1957;  Irwin,  1957;  Kanninen  and  Popelar,  1985) 

a.j(r,e)=  I,j(0)-KTij(e)-^ 

where  (0)  functions  are  determined  from  eqs  (2. 2), (2. 3).  To  specify 
functions  (9)  and  H'-j(9)  ,  consider  formulation  of  the  plane  problem  in 

terms  of  Airy's  stress  function  of  complex  variable  z  =  x  +  iy  =  re  . 
Following  to  Kanninen  and  Popelar  (1985),  one  can  write 

+  Oyy  =  2ia’(z)+  n'(5)l  , 

Oyy  ~  -  2i  =  2  I  z  n"  (^  +  m"  (z)  1  , 


r  A'P..(0)-i-...  (2-16) 

Vh  ^ 


(2.17) 


where  £2(z)  and  co(z)  are  holomorphic  functions,  the  overbar  denotes  the 
complex  conjugate.  For  Mode  II,  the  stress  function  is  an  odd  one  of  the  y- 
coordinate,  and  £2  and  co’  take  the  form  (Kanninen  and  Popelar,l985) 

Q(z)  =  iaz^^'  ,  03'(z)  =  ibz^"^^  (2.18) 

where  a,b  and  X  are  real  constants  to  be  determined  from  the  boundary 
conditions;  X>-1. 

After  substitution  of  eqs  (2.18)  into  eqs  (2.17),  one  can  find 

^xx  ^yy  -  -4a(X+l)  sin  XB  ,  (2.19) 

-  itj^  =  -(X+l)  r^  {  (2a+b)sin  A0+  aX  sin  (X-2)6 

+  i  [  aX.  cos  (A-2)0  +  b  cos  A0  1  }  (2.20) 

For  traction-free  crack  surfaces,  eq.(2.20)  yields  the  boundary  conditions 

(2a  +  b  +  aA)  sin  Arc  =  0  , 

(aX  +  b)  cos  Arc  =0  ,  (2.21) 

which  result  in  eigenvalue  solutions  of  the  boundary  value  problem  as 

sin  Arc  =  0  or  cos  Arc  =  0  , 

or,  equivalently, 

X  =  n/2  with  n  =  -1,  0,  1,  2,  ... 

For  the  first  eigenvalue,  X=-l/2,  eqs  (2.19)-(2.21)  lead  to  the 
singular  stress  field  given  by  eqs  (2.2).  The  second  eigenvalue,  A=:0, 
corresponds  to  the  case  in  which  all  stress  components  vanish,  i.e. 

=  Tyy  =7,y  =  0  (2.22) 

Thus,  in  contrast  to  Mode  I,  the  crack  tip  stress  field  in  Mode  II  does  not 
contain  constant  value  terras  (T-terms). 

Substituting  A=l/2  into  eqs  (2.19)-(2.21),  one  can  find  'P-terms: 

=  sin- (2  +cos^-)  , 

2  2 

=-sin^cos^y  ,  (2.23) 

'J'xy  =cos|  (1  • 

or,  in  polar  coordinates, 

'P_.  =  sin—  (2  +  5  cos—  cos-— )  , 

”  2  2  2 

=  -5sin—  cos^  .  (2.23a) 

ui  5  -  30  3  e 

'P.fl  =  -COS0COS - cos— 

“2  2  2  2 

For  certain  geometry  and  loading  conditions  the  amplitude  A  of  non¬ 
singular  4* -components  of  the  stress  field  (2.16)  can  be  determined  from 
results  of  numerical  treatment  of  a  boundary  value  problem.  For  example, 


by  analogy  with  approach  of  Larsson  and  Carlsson  (1973),  one  can  w-rite 

A  =  UTn^_^Q  aJ  —  I  (r,e=0)  -  (r,e=0)  ] 

y  ^  ^appi 

where  and  are  stress  components  for  the  actual  geometry  of  a 

cracked  body  and  for  geometry  corresponding  to  the  boundary  layer 
formulation  with  one-parameter  (Kjj)  characterization  of  the  crack  tip 
fields. 

An  approximate  assessment  of  the  'V  -stress  amplitude  can  be 
obtained  for  particular  crack  geometries  as  follows.  Consider  a  plate  of 
width  W  containing  an  edge  crack  of  length  u;  let  H=W-a  be  a  ligament  size. 
The  plate  is  subjected  to  shear  stresses  tappi  acting  parallel  to  the  crack 
line.  With  terms  of  order  r'  be  omitted  in  eq.(2.16),  the  equilibrium 
condition  yields 

H 

(r,e=0)  dQ  =  r^ppi  W  , 

0 

and.  after  integrating, 

A  =  l(  W  _  JiL^/T) 

2  H  T,pp,  V 

Making  use  of  conventional  representation  of  the  stress  intensity  factor  in 
the  form 

“^appi  ^ 

where  Tjj  is  a  geometrical  factor,  we  get  finally 


-|( 


- ) 
H 


(2.24) 


Consider  two  examples  of  calculating  the  non-singular  stress  term 
amplitude  for  different  crack  geometries  in  Mode  II. 


(1)  Edge  crack  of  length  a-H=W/2  in  a  plate. 

-According  to  results  published  by  Tada  et  al.  (1973), 

Kn  =  1.252  , 

and  formula  (2.24)  yields 

A  =  0.344 

(2)  Semi-infinite  crack  in  an  infinite  plate. 

This  example  is  a  limiting  case  of  the  previous  one  with  .  The 

geometrical  factor  equals  (Savruk,  1988) 

Y 

rjj  - -  , 

and  one  can  find  from  eq.(2.24) 

A  =  1.09 
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For  the  latter  geometry  an  analysis  of  the  plastic  zone  evolution  with 
variation  of  the  stress  level  is  performed  below.  The  Huber-Mises  yield 
condition  (2.12)  can  be  rewriten  in  the  form 

I  [(<yxx -^yy^  +  ,  (2.12a) 

or,  equivalently,  for  plane  strain  conditions 

(I  -  V  +  ~  'i(a^Cyy-x^^)  =  3x^  (2.^6) 

As  it  is  mentioned  above,  the  yield  condition  for  plane  stress  may  be 
derived  after  that  for  plane  strain,  eq.(2.12b),  putting  v=0. 

Introducing  dimensionless  parameters 

a^j  =  Ojj  /  Xy  ,  R  =  r(XY  /K^)  ^  , 

the  stress  components  and  the  yield  criterion  are  of  the  form 

aij(R,e)  =  — i=-  Sij(0)  +  AKn  (^)2  'I'ijO)  ,  (2.25) 

V  27cR  V 

(1  -  V  +  v^  )  (a;f;c  +  Qyy  f  -  3  (oxx  ^ yy  *  )  =  3  (2.26) 

Equations  (2.25)  and  (2.26)  predict  dependence  of  the  dimensionless  radius 
of  the  plastic  zone, 

Rp  —  tp(XY  /Kjj)  ^  , 

on  the  applied  stress  level;  this  conclusion  is  in  agreement  with  the  results 
of  Larsson  and  Carlsson  (1973)  and  Rice  (1974)  on  Mode  I. 

The  closed  form  solution  of  eqs  (2.25)  and  (2.26)  for  Rp  is  rather 
complicated,  variation  of  Rp(0)  has  been  derived  by  numerical  treatment  of 
eqs  (2.25)  and  (2.26)  for  different  stress  levels  and  two  values  of  Poisson's 
ratio:  v=:0  and  v=0.33.  Results  for  both  plane  stress  (v=0),  Fig.5a,  and  plane 
strain  (v=0.33).  Fig. 5b,  indicate  significant  change  in  the  shape  of  plastic 
zone  as  compared  to  shapes  derived  from  the  model  with  one-parameter 
characterization  of  the  stress  field  near  a  crack  tip  (Fig.4).  The  present  two- 
parameter  solution  predicts  more  rapid  growth  of  the  plastic  zone  ahead  of 
the  crack  tip  than  behind  it.  With  increasing  the  load  level,  the  major 
plastic  zone  area  concentrates  ahead  of  the  crack  tip,  which  is  more 
characteristic  of  the  plane  strain  case.  For  the  latter  the  zone  hight 

ahead  of  the  crack  tip  is  larger  that  behind.  These  features  of  rhe  plastic 
zone  behavior  at  the  crack  tip  in  Mode  II  differ  from  those  predicted  by 
the  one-parameter  model  (Fig.4)  and  agree  with  the  finite-element  results 
by  Banks-Sills  and  Sherman  (1990). 

Note  that  previous  solutions  for  the  plastic  zone  geometry  do  not 
incorporate  Irwin's  correction  for  plasticity.  Since  conditions  s=  (Syy  =  0, 
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0  are  hold  at  9=0,  Irwin's  approach  is  directly  applied  to  evaluation 
of  the  plastic  zone  size  ahead  of  the  crack  tip.  This  results  in  the  value  of 
rp(0=O)  twice  larger  than  predicted  before  (Figs  4,5),  while  the  values  of 
rp(6=7t)  and  the  plastic  zone  hight  do  not  change.  It  is  briefly  discussed  in 
Chapter  6  which  assessments  of  the  plastic  zone  size  and  shape  are  more 
realistic. 

2.6.  Stress  Intensity  Factors  for  Cracks  Subjected  to  In-Plane 
Shear  Loading 

Within  the  SSY  approximation,  the  size  and  shape  of  plastic  zone  near  a 
crack  tip  are  expressed  via  the  stress  intensity  factor.  Some  reference  K^- 
solutions  for  plane  crack  problems  are  summarized  in  this  section.  Note 
that  in  stress  intensity  factor  analyses  actual  boundary  conditions  are  often 
formulated,  according  to  the  superposition  principle  (Bueckner,  1958),  in 
terms  of  self-equilibrium  stresses  applied  to  crack  surfaces.  With  such 
transformation  the  formulation  of  a  boundary  value  problem  is  simplified 
and  some  efficient  techniques  for  calculation  of  the  stress  intensity  factors 
can  be  employed  (e.g.,  Green's,  or  weight  function,  method). 

(1)  Central  crack  of  length  2^?  in  an  infinite  plate  (Fig.6.1). 

For  arbitrary  stress  distribution  on  crack  surfaces,  t(£),  stress  intensity 
factors  at  right  (+)  and  left  (-)  crack  tips  are  determined  as  follows  (Savruk, 
1988) 

a 

Kjf  =-p=-  J  x(^)-\I  (a  +^)/(a  -  ^)  d^  , 

ytta  _a 

a 

=-^  J  x{%)^lia-ty{a  +  ^)  d^ 

'yixa 

In  the  case  of  constant  stress  distribution,  x(^)  =  Xq  “  const,  eqs  (2.27)  yield 

iC^  =  /Cn=  XoVjt7  (2.28) 

If  a  pair  of  forces  ±Q  acts  at  points  (b,0*)  and  (b,0~)  respectively,  the  stress 
intensity  factors  are 

K^  =  -^'^ia  +  b)/ia-b)  ,  ~  (a-b)/(a  +  b)  (2.29) 

V  Tca  y  Tca 

For  constant  stresses  Xq  ,  applied  near  the  crack  tip  regions,  bi\^  jia  (as 
in  Dugdale's  model),  one  can  find 

Kn-  ^-XQ^fm  (-  -  sin-1 

It  2  a 


(2.30) 


Cracks  of  length  2a  with  a  distance  d  between  cracks  are  considered.  The 
Kn'Solution  for  uniform  stress  distribution,  Xq  ,  on  crack  surfaces  (or 
uniform  remote  stresses  Xq  )  has  been  derived  by  Savruk  (1988): 

iCn  =  ^K(th^)  A  /  (a/3>.)  th(^-^  (2.31) 

jt  4  >/  4 

with  X-2ald,  and  K  is  the  complete  elliptical  integral  of  the  first  kind.  The 
solution  for  an  arbitrary  distribution  of  shear  stresses  is  rather 
complicated  and  may  be  found  in  Savruk’ s  (1988)  book, 

(3)  Semi-infinite  crack  in  an  infinite  plate  (Fig. 6. 31. 


Stress  intensity  factor  for  an  arbitrary  stress  distribution  is  determined  as 
follows  (Panasyuk  et  al.,  1976): 


(2.32) 


The  function  x(^)  is  assumed  in  such  a  way  that  the  integral  in  eq.(2.32) 
takes  a  finite  value.  If  the  region  near  crack  tip,  0  i  ^  s.  a,  is  subjected  to 
constant  shear  stresses  Xq,  eq,(2.32)  yields 


2V2 

7t 


(2.33) 


In  the  case  of  a  pair  of  shear  forces  +Q  acting  at  points  (fc,0+)  and  (6,0  )  one 
can  find 


ifn  =  Q 


(2.34) 


(4)  Periodic  array  of  semi-infinite  cracks  in  an  infinite  plate  (Fie. 6.4). 


An  approximate  solution  of  this  problem  has  been  derived  in  a  closed  form 
by  special  approximation  of  a  kernel  of  the  singular  integral  equation 
(Savruk,  1988) 


t(^) 

V 1  -  exp(-3x^/d) 


0 


(2.35) 
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(5)  Edge  crack  of  Finite  length  a  in  a  half-planei  (Fig, 6. 51 


Several  approximate  solutions  of  this  problem  have  been  derived  in  the 
literature.  According  to  Panasyuk  ei  al.  (1976),  the  stress  intensity  factor 
for  an  arbitrary  stress  distribution  is  of  the  form 


_  llaa^  J°  x(^)  d^ 

V  ^  0 


(2.36) 


with  a=2xV(i=^-4).  For  the  stress  distribution  of  polynomial  type,  t(^)  = 
,  eq.(2.36)  is  reduced  to 


^ ^ —  (2.37) 

°  V  a  r(n/2+l) 

Making  use  of  the  alternating  technique,  Hartranft  and  Sih  (1973) 
have  found  Green’s  function  which  is  the  value  of  the  stress  intensity 
factor  due  to  a  pair  of  shear  forces  ±Q  acting  at  points  (c,0'^)  and  (c,0‘)  of 
crack  surfaces  respectively 

kn  =  [  1  +  (1  -  e^)  (0.2945  -  0.3912e^  +  0.7685e^ 

-  0.9942e®  +  0.5094e®)  ]  /  (2.38) 


with  E-c/a.  The  stress  intensity  factor  of  an  arbitrarily  loaded  edge  crack  in 
a  half-plane  can  be  calculated  by  integrating  Green's  function  (2.38)  over 
the  crack  length.  For  example,  for  constant  shear  stresses  one  can  get 

Kn  =  1.1215  (2.39) 

At  E=0,  eq.(2.38)  yields 

iCn  =  1.2945  ^  (2.40) 

"Vita 

which  is  the  solution  of  the  problem  when  two  opposite  shear  forces  Q  are 
applied  at  the  crack  mouth  (^=0,  ti=±0). 

If  a  tensile  load  P  parallel  to  the  ^-axis  is  applied  to  a  half-plane 
boundary  in  the  point  (0,b)  (Fig. 6.5),  the  corresponding  value  of  the  stress 
intensity  factor  was  derived  by  integrating  (2.38)  (Rooke  and  Jones,  1979) 

Ka  -  -  P  (1.294  -  I.184(;  -H  5.442J;^  -  28.14!;^  -h  41.80C'‘ 

YJtn 

-  22.38i;5 -t- 3.162C*)  (2.41) 


with  l^-bl(b+a). 

Freund  (1978)  has  solved  the  problem  with  prescribed  displacements 
of  a  half-plane  boundary 

u^  =  5  =  const  for  ^=0,  ti>0  , 

=  -6  =  const  for  ^=0,  •n<0  ; 


the  result  is 
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/Cn=:-  (2.42) 

where  jj.  is  the  shear  modulus,  ic=:3-4v  for  plane  strain  and  k=(3-v)/(1+v)  for 
plane  stress  conditions. 

(6)  Edge  crack  in  a  half-plane  with  a  stiffened  edge  under  uniform 
tension  (Fig.6.6). 

Let  a  rigid  stamp  of  length  b  be  joined  with  a  half-plane  edge  as  shown  in 
Fig.6.6.  The  force  P  parallel  to  the  crack  line  acts  on  the  half-plane  through 
the  stamp.  This  problem  has  been  solved  by  Hasebe  (1979,  1981)  by 
conformal  mapping  technique.  Two  types  of  boundary  conditions  were 
studied;  (a)  loading  with  no  stamp  rotation  and  (b)  loading  with  stamp 
rotation  at  an  angle  cp;  corresponding  results  are  presented  in  the  table 
below  (ifj  -values  are  given  over  the  lines,  -  values  are  given  below  the 
lines). 


aJb 


Kij^/iP^I  Ti/a) 


/Cjjj  /(pcp-s/jcT) 


K 


K 


3.0 

2.5 

2.0 

5/3 

1.0 

3.0 

2.5 

2.0 

5/3 

1.0 

0.002 

0.020 

0.019 

0.017 

0.015 

0.W7 

0.004 

0.005 

0.006 

0.006 

0.006 

0.001 

0,681 

0.709 

0.715 

0.692 

0,454 

0.178 

0.231 

0.294 

0.341 

0.413 

0.1 

0.0S9 

0J)84 

0.077 

0.071 

0.053 

0.238 

0.258 

0.285 

0.308 

0.379 

0.049 

0.050 

0.049 

0  049 

0.044 

0.228 

0.253 

0.285 

0.313 

0.392 

0.4 

0.152 

0.150 

0.146 

0.143 

0.134 

0.097 

0.110 

Q.128 

Q.145 

0.202 

0.115 

0.114 

0.113 

O.lll 

0.105 

0.158 

0.178 

0.205 

0.229 

0.307 

0.8 

0,207 

0.206 

0.205 

0.205 

0.203 

0.038 

0.043 

0.050 

0.057 

0.080 

0.184 

0.184 

0.183 

0.183 

0.181 

0.100 

0.113 

0.131 

0.147 

0.200 

1.0 

0,224 

0.224 

0.223 

0.223 

0.223 

0024 

0.027 

0.032 

0.036 

0.051 

0.212 

0.213 

0.213 

0.213 

0.212 

0.080 

0.090 

0.105 

0.118 

0.161 

5/3 

0.251 

0.251 

0.251 

0.251 

0.252 

0.005 

0.006 

0.006 

0.007 

0.010 

0.279 

0.280 

0.281 

0.282 

0.284 

0.033 

0.037 

0.043 

0.048 

0.065 

2.5 

0.200 

0.260 

0.260 

0.260 

0.261 

0,001 

0.001 

0.001 

O.OQl 

0.002 

0.323 

0.324 

0.325 

0.325 

0.327 

0.015 

0.016 

0.019 

0.021 

0.029 

10. 

0.262 

0.263 

_D.263 

0.263 

0.263 

0. 

0. 

0. 

0. 

JL 

0.390 

0.390 

0.391 

0.391 

0.392 

0.001 

0.001 

0.001 

0.001 

0.002 

In  both  cases  the  mixed  Mode  I-II  crack  tip  deformation  takes  place.  When 
the  stamp  and  the  crack  are  of  equal  length,  a!b  =  \,  in  case  (a)  the  absolute 
values  of  Kj  andKjj  are  close,  while  in  case  (b)  the  ratio  1/Cjj/iCjl  is  of  the 
value  about  3. 

(7)  Periodic  array  of  parallel  edge  cracks  in  a  half-plane  (Fig-6.7'). 

Numerical  results  of  solving  this  problem  by  the  method  of  singular 
integral  equations  (Panasyuk  et  al..  1976)  are  summarized  below  for  the 


case  of  constant 

shear 

stresses 

Xq  applied  to  crack 

surfaces: 

aid 

0. 

0.2 

0.4  0.6  0.8 

1.0  2.0 

3.0 

iCn/ftQ  V  7t<3  ) 

1.1214 

1.1320 

1.2072  1.3291  1.4575 

1.5797  2.0941 

2.5075 

(8)  Laver  of  finite  length  bounded  to  a  half-space  (Fig. 6.8). 

This  problem  has  been  solved  by  the  method  of  singular  integral  equations 
(Chantaramungkorn  and  Keer,  1975).  Both  the  layer  and  half-space  are 
assumed  to  be  of  identical  material.  Constant  pressure,  p,  acts  parallel  to 
the  crack  plane.  .Numerical  /Cj  and  values  are  presented  in  the  table 
below  for  two  crack  tips  (A  and  B): 


c_ 

h 

d 

c 

pN  TtC 

pv  zc 

fC,(A) 

p\  ;tc 

P'i  TIC 

0.5 

1.25 

0,648 

0.458 

1.349 

-1.309 

l.O 

1.25 

0.367 

0.188 

0.398 

-0.325 

1.50 

0.385 

0.198 

0.390 

-0.321 

2.00 

0.387 

0.200 

0.390 

-0.320 

3.00 

r ' 

0.200 

0,390 

-0.320 

2.0 

1.25 

C 

0.066 

0.188 

-0.076 

3.0 

1.25 

0.185 

0.033 

0.145 

-0.033 

Note  that  practically  for  all  values  of  geometrical  parameters  studied,  the 
condition  i/CQ/iCjj  <  1  holds. 
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(9)  Edge  crack  in  a  strip  (Figi.6.9). 

Approximate  formula  for  the  stress  intensity  factor  at  constant  shear 
stresses,  Xq,  applied  to  crack  surfaces  is  given  by  Tada  et  al.  (1973): 

^n  =  '^oV^/o(“)  . 

/q(^)  =  [  1.122  -  0.561(^  +  0.085(^2  +  0.1 80(:^)^]  /  ^ ^  (2.43) 


which  is  said  to  be  accurate  within  2%  for  0  <  ^  <  1 .  In  the  case  of  a  pair  of 

shear  forces  Q  acting  at  the  crack  edge  (see  Fig. 6. 9)  an  expression  for  is 

of  the  form  (Tada  et  al.,  1973) 

An  =  t  1  -30  -  0.65(^)  +  0.37(^)  2+  0.28(^)  ^  ^  1  -  (2.44) 

which  is  accurate  within  1  %  for  0  <  ^<  1 . 

W 

Making  use  of  an  asymptotic  interpolation  technique.  Cheng  and 
Finnic  (1990)  have  derived  the  stress  intensity  factor  for  an  arbitrary 

stress  distribution,  x(^)  =  Tq  t(^),  as  follows: 


Kn  =  Tq  >/it^/o(:|r)  . 


W' 

with  /o(^)  to  be  determined  from  eq,(2.43),  and 
a/W 


(2.45) 


=  Jal['  -  ■ 


(2.46) 


(10)  Edge  crack  in  a  three-ooint  bend  specimen  (Fig.6.1Q). 

Mixed-mode  stress  intensity  factors,  and  Kjj,  for  an  asymmetrically 
located  edge  crack  in  a  three-point  bend  specimen  have  been  calculated  by 
the  boundary  collocation  method  (Ke  et  al.,  1978).  Dimensionless  functions 
Fj  = /Cj /[Af/(rW  ^^)]  and  =  Kjyl[QI{tW  are  summarized  in  the  table 

below  with  Fj  being  above  the  line  and  Fq  being  below  the  line  (Af  and  Q 
are  the  bending  moment  and  shear  force  in  a  cracked  section;  r  and  W  are 
the  specimen  thickness  and  width); 


a 

W 

2BJL 

0 

1/6 

2/6 

3/6 

4/6 

5/6 

11/12 

0.40 

7.71 

R..S0 

g55 

8.36  . 

8.33 

SJ£L  . 

0 

1.032 

1.400 

1.350 

1.298 

1.376 

1.644 

0.45 

^.S6 

9  "7 

9.38 

948 

*  »* 

0 

1.142 

1.562 

1.488 

1.466 

1.464 

... 

0.50 

10.27 

11. 4S 

11. -0 

11-60 

ILU_ 

11.59 

U.5.3 

0 

1.410 

i.S64 

1.840 

1.664 

1.660 

1.760 

0.55 

12.11 

13.30 

Ii6(L 

13.03 

12.90 

13.46 

•  •  • 

0 

1.588 

1.980 

2.050 

1.976 

2.100 

... 

0.60 

14.47 

14,25 

14^ 

14  91 

\±M 

14.74 

14,50 

0 

2.348 

2.248 

,  2.276 

2.320 

2.294 

2.090 

The  value  of  /Cj  is  found  to  be  much  higher  than  iCg  . 

(11)  Edge  crack  in  a  four-point  bend  specimen  (Fig. 6.11). 

Pure  Mode  II  crack  tip  deformation  is  expected  in  the  loading  scheme 
presented  in  Fig. 6. 11.  This  specimen  has  been  employed  by  Tohgo  and  Ishii 
(1992)  for  measurements  of  the  fracture  toughness,  .  The  stress 

intensity  factor  has  been  found  from  finite-element  calculations: 

iCj  =  1.37  ,  F  =  -P  (2.47) 

“  Wt  8 

where  W  and  f  are  the  specimen  width  and  thickness,  respectively.  Note 
[hat  varying  distances  in  the  loading  frame,  combined  Mod?  l-II  conditions 
can  be  obtained  with  various  mixity  factor  (Tohgo  and  ishii,  1992). 

(12)  Edge  crack  in  a  cantilever  beam  (Fig,6.12). 

Stress  intensity  factors,  Kj  and  Kq  ,  have  been  found  by  the  body  force 
method  (Murakami,  1980),  Dimensionless  functions 

Fi  =  Kj  /[(6rL/lV2)V^]  .  Fji  =  /Cn  /[iT/W)i 
(7’  =  'qH0  are  presented  below  for  different  values  of  the  crack  depth: 


2  4 


alW 

0.1 

0.2 

0.3 

0.4 

0.5 

Fi 

0.94 

1.00 

1.12 

1.25 

1.36 

Fq 

0.43 

0.68 

0.94 

1.17 

1.36 

Since  LI'W>\,  then 


\ 

w 


»  1, 


i.e.  Mode  I  deformation  is  predominant  for  this  specimen. 

(13)  Edge  crack  in  a  double  cantilever  beam  (Fig.6.13). 

A  double  cantilever  beam  is  considered  with  different  cross  sections, 
i4i  =  t/iiand  A2=th2  ,  where  r  is  the  beam  thickness  (Fig. 6. 13).  A  pair  of 
forces.  P.  is  applied  to  the  beam  edges.  The  stress  intensity  factor  is 
determined  as  follows  (Tada,  1974): 

^  (2.48) 


V2rA 


with 


A  ~ 


C  =  1 
C  = 


Vl-v' 


Aj  +A2 

for  plane  stress, 
for  plane  strain. 


If  Ai=A2  and  a  horizontal  displacement  5  of  the  beam  is  given,  eq.(2.48)  is 
reduced  to  the  following: 

£5  /2J 


a 


t 


(2.48a) 


(14)  Ead-nowhgd  flgjmrs  ?p6£imfin.,in-thrc.6--pQini  .hsndiitf 

This  type  of  specimen  has  been  employed  by  Mall  and  Mol  (1991)  in  Mode 
II  fracture  toughness  testing  of  ceramic  composites.  In  preparation  phase, 
the  specimen  is  notched  on  one  end  and  subjected  to  Mode  I  loading  to 
create  a  sharp  crack  of  length  a  which  is  shown  in  Fig.6.14.  The  fracture 
toughness  is  expressed  in  terms  of  the  critical  energy  release  rate  as 
follows: 


G  ~^C  dC 
2f  da 


(2.49) 


where  Pc  is  a  critical  force,  t  is  the  specimen  thickness.  C  is  the  specimen 
compliance.  An  approximate  formula  for  evaluation  of  the  C  value  is  given 
by  Russel  and  Street  (1982): 


(2.50) 

iEl(.h/L)3 

However,  Mall  and  Mol  (1991)  demonstrated  that  the  real  comliance 
function  differs  essentially  from  that  determined  by  eq.(2.50).  and  thus  C- 
calibration  curves  are  to  be  found  experimentally  for  certain  specimen 
geometry  and  test  conditions  (temperature  etc.). 


(15)  Compact  specimen  in  shearing  (Fig.6.15) 


The  stress  analysis  of  this  specimen  has  been  performed  by  the  finite- 
element  method  (Hoyniak  and  Convvay,  1979).  Two  types  of  boundary 
conditions  (with  and  without  restrictions  on  horizontal  displacements  of 
the  applied  load)  and  two  different  directions  of  the  applied  load  have 
been  considered;  these  are  refered  to  as  Cases  1  to  4  in  Fig.6.15.  In  general, 
mixed  mode  type  of  the  crack  tip  deformation  is  observed  for  this 
specimen.  Figure  7a  shows  variation  of  the  stress  intensity  factor  for  the 
Case  2  of  loading  found  by  the  finite-element  method  (Hoyniak  and 
Conway.  1979)  and  the  boundary  collocation  method  (Jones  and  Chisholm. 
1975),  for  three  values  of  the  characteristic  size  of  the  specimen:  H  =  0.5, 
l.O  and  1.5  in.  In  Fig. 7b.  both  the  Mode  II  stress  intensity  factor,  ,  and 

the  mixity  factor,  5  =  |/Cjj  |  /  "s/  .  are  presented  for  different 

schemes  of  loading  (after  Savruk.  1988).  One  can  conclude  that  nearly  pure 
Mode  II  crack  tip  deformation  is  found  for  the  Case  3. 

(16)  Punch-through  shear  specimen  (Fig  6,16). 

This  specimen  (cubic  in  shape)  has  been  applied  for  fracture  toughness 
(Kj^q)  measurements  of  concrete,  i.e.  extremely  brittle  material  (Davis  et  al., 
1985).  The  calibration  of  the  stress  intensity  factor  has  been  performed  by 
the  finite-element  method.  The  specimen  dimensions  in  that  analysis  were: 
W  =  100  mm,  b  =  10  mm.  H  =  30.  40  and  50  mm.  Numerical  results  by 
Davis  et  al.  (1985)  on  the  dimensionless  stress  intensity  factor 


l'ii  = 


2W{W-a-b) 

P 


^Ul 


for  a  larger  crack  (of  length  a)  are  presented  in  the  table  below. 
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aAV 

for  H=50  mm 

Fu  for  H=40  ram 

Fji  for  H=30  mm 

0.25 

4.833 

4.151 

3.589 

0.30 

4.802 

4.109 

3.573 

0.35 

4.769 

4.088 

3.550 

0,40 

4.743 

4.071 

3.521 

0.45 

4.690 

4.067 

3.510 

0.50 

4.394 

3.809 

3.335 

0.55 

4.438 

3.863 

3.405 

0.60 

4.544 

3.978 

3.484 

For  prescribed  specimen  dimensions  the  opening  mode  stress  intensity 
factor  was  found  to  be  of  the  order  lower  than  the  value  of  Kq. 


(17)  Compact  Mode  II  fracture  specimen  (Fig.6.17). 

This  specimen  has  been  designed  and  employed  for  fracture  toughness 
(Kqc  and  measurements  of  perspex  (extremely  brittle  material) 

(Banks-Sills  and  Arcan,  1986)  and  A1  7075-T7351  (typically  ductile 
material)  (Banks-Sills  and  Sherman,  1991).  The  specimen  geometry  is 
presented  in  Fig. 6. 17.  The  load  is  applied  through  a  special  frame  so  that 
no  bending  moment  acts  at  a  cracked  section  (this  is  shown  schematically 
in  Fig.6.17).  The  finite  element  method  has  been  used  for  elastic  and 
elastic-plastic  analysis  of  the  specimen.  An  approximate  formula  for  the 
stress-intensity  factor  was  derived  by  Banks-Sills  and  Arcan  (1986).  Also 
calibration  curves  for  evaluation  of  the  Jq-vuIuc  are  presented  by  Banks- 
Sills  and  Sherman  (1991). 

It  should  be  noted  in  conclusion  of  this  brief  review  on  Mode  II 
stress  intensity  factors,  that  application  of  cracked  geometries  shown  in 
Figs  6.11,  6.14,  6.15  (Case  3),  6.16  and  6.17  allows  for  pure  Mode  11  crack 

tip  deformation.  For  the  other  geometries  mixed  Mode  I-II  conditions  take 

place  with  different  ratio  of  Kj  to  Kjj  (depending  on  a  geometry  and 

scheme  of  loading).  It  is  important  to  note  than  at  iCj  <0  compressive 

stresses  act  along  a  crack  line  inducing  contact  and  friction  between  crack 
surfaces;  this  results  in  reducing  an  actual  in-shear  mode  stress  intensity 
factor  and  leads  to  necessity  to  reformulate  boundary  conditions.  Also,  in 
the  presence  of  the  crack  surfaces  interaction,  the  proof  of  the  7-integral 
path-independence  becomes  invalid  and  application  of  the  /-integral  to 
ductile  fracture  analysis  is  questionable. 
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2.7.  Some  Results  on  Mode  II  Fracture  Toughness  Measurement 

Some  results  on  determination  of  fracture  toughness  in  Mode  II  available 
from  the  literature  are  summarized  in  Table  2.1.  For  comparison  purposes, 
values  of  Mode  I  fracture  toughness,  KhqIKiq  ratio  and  crack  initiation 
angle  0j(with  respect  to  the  original  crack  line)  are  listed. 

Table  2.1 


Reference  Material 


MPaV' m 


Banks-Sills  and 

Arcan  (1986) 

perspex 

Banks-Sills  and 
Sherman  (1991) 

A17075-T7351 

29.2 

Jones  and 

Chisholm  (1979) 

Al  2024-'^4 

37.4 

Hoyniak  and 
Conway  (1979) 

Al  2024-T4 

37.4 

Tohgo  et  al.  (1990) 

steel  SM41A 

500’ 

Tohgo  and  Ishii 
(1992) 

.Al  6061 -T  651 

15.3^ 

Davies  et  al.  (1985) 

soil  cement 

Awaji  and  Sato 
(1978) 

graphite  7474 
graphite  SMI -24 

0.94 

0.81 

Mall  and  Mol 
(1991) 

alumosilicate 

ccratmc 

composite 

Suresh  ec  al. 
(1990) 

AliOj  ceramic 

3.32: 

8j  Specimen 

ATjc 


MPaVm 

(deg) 

0.89 

63-70 

Fig.  6.17 

54.1 

1,85 

0 

Fig.  6.17 

44.0 

1.18 

=0 

Fig.  6,15 

55.8 

1.49 

=0 

Fig.  6.15 

205’ 

0.64 

0 

Fig.  6.11 

45.2t 

1.72 

0 

Fig.  6.11 

0.42 

=0 

Fig.6.16 

1.09 
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3.  Plastic  Stress  and  Strain  Fields  near  a  Crack  Tip 
in  Mode  II 

In  previous  chapter,  analyses  of  the  crack  tip  plastic  zone  size  and  shape 
are  based  upon  the  linear-elastic  solution  for  the  stress  field.  The  only 
plasticity  parameter  involved  in  those  analyses  is  the  yield  stress,  Xy-  In 
reality,  there  exists  material  strain  hardening  which  affects  significantly 
material  behavior  in  yielding  and  leads  to  redistribution  of  stress  and 
strain  fields. 

The  dominant  singularity  at  a  crack  tip  in  a  power  strain  hardening 
material  has  been  investigated  by  Hutchinson  (1968a)  and  Rice  and 
Rosengren  (1968)  (HRR-solution).  In  those  papers,  Mode  I  stress-strain 
fields  in  plane  stress  and  plane  strain  have  been  derived  for  material  of 
the  Ramberg-Osgood  type,  as  well  as  for  perfectly  plastic  one.  Further 
studies  were  devoted  to  the  elastic-plastic  analysis  of  the  crack  tip  fields  in 
pure  Mode  II  plane  strain  (Hutchinson,  1968b;  Shih,  1974)  and  pure  Mode 
III  (Rice,  1968b),  as  well  as  combined  Mode  I-II  (Shih,  1974;  Nemat- 
Nasser  and  Obata,  1984),  I-III  (Pan  and  Shih,  1989)  and  II-III  (Pan  and 
Shih,  1990). 

Note  that  finite-element  analyses  of  mixed-mode  crack  tip  fields  led 
to  some  important  conclusions:  (1)  the  dominant  HRR-singularity  derived 
originally  within  deformation  plasticity  is  valid  when  incremental 
plasticity  is  employed  to  describe  material  behavior;  (2)  in  mixed-mode 
conditions,  there  is  only  slight  deviation  of  the  dominant  stress  field 
singularity  order  from  that  in  the  HRR-solution  (Pan  and  Shih,  1990).  Thus, 
the  HRR-solution  appeares  to  be  a  good  approximation  of  real  stress  and 
strain  fields  near  a  crack  tip  in  metallic  materials. 

In  what  follows,  a  brief  review  of  the  HRR-solution  for  pure  Mode  I 
and  Mode  II  is  presented.  In  addition,  variation  of  the  dominant  stress  and 
strain  fields  at  a  crack  tip  in  Mode  II  versus  strain  hardening  exponent  is 
analysed  both  for  plane  strain  and  plane  stress. 

3.1.  Dominant  Singularity  of  the  Crack  Tip  Fields 

A  plane  crack  is  considered  in  an  elastic-plastic  material  obeying  the  J2- 
flow  theory.  The  stress-strain  curve  in  shearing  is  approximated  by  the 
Ramberg-Osgood  equation 

T  =  (3.1) 

where  Xy  and  Yy='^y/^  yield  stress  and  strain  respectively,  p.  is  the 

shear  modulus,  N  is  the  strain  hardening  exponent  (N  >1),  a  is  the  material 
constant.  The  plastic  deformation  is  assumed  to  be  incompressible, 
independent  of  the  hydrostatic  stress  component,  a^^/3,  and,  thus,  be  a 
function  of  the  deviatoric  stress  components,  sjj,  i.e. 
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£. .  =  e  ®  +  £  p  , 

ij  ij  u  ’ 


(3.2) 


where 


e  1+v  ,  l-2v  5 

E  ^  3E 


e.P  =  s  ■ 

2  a  ’ 


(3.3) 


with 


^ij  -  2  ^kk^ij 


(3.4) 


The  effective  stress,  ,  and  plastic  strain,  eP  ,  corresponding  to  the  J2-flow 
theory  are 

a,  =  VilT 


^.=  •^75^^  Sij 


or,  in  explicite  form. 


=  I  “[(ctxx-<Jyy)^+(<^yy-^z2)^+(<^22-Ctxx)‘]  +  3(x^^+ty2+x^2)|  1/2  ^  (3  5^) 


eP  =  A  /-£  p£  P 

«  V  ;  iJ  ij 


=  [(=,!)'+  (Sy")'*  2(e,f)'+  2(e„P)=]  (3.5b) 


Making  use  of  eqs  (3.1)-(3.3),(3.5)  the  following  relation  is  obtained 
between  the  effective  stress  and  plastic  strain: 


eP  =  a£Y 


(3.6) 


with  eY=7Y/V^.  cJy=V3  Ty.  Since  the  near  crack  tip  fields  with  a  high  level  of 
plastic  deformations  are  considered,  the  stress-strain  relation  may  be 
approximated  by 

e-E,j>  =  |ae,[i')''-‘^  (3.7) 

2  I  Oy  J  Oy 


where  eqs  (3.3)  and  (3.6)  are  taken  into  account. 

The  order  of  the  dominant  stress  and  strain  fields  singularity  can  be 
obtained  directly  from  the  analysis  of  the  path-independent  /-integral 


( 


30 


(Rice,  1968a).  Thus,  for  a  circular  contour  of  radius  r  enclosing  the  crack 
tip  the  J-integral  (2.5)  is  expressed  as  (Rice  and  Rosengren,  1968) 


r  Piti. 

J  =  r  J  (W  cosd  -  CijUj  — ^  )  d0 


(3.8) 


-7t 


Since  the  value  of  J-integral  is  independent  on  r,  both  terms  in  brackets  of 
eq.(3.8)  have  to  reveal  the  singularity  of  the  form 

-1  _  dui  (3  9) 


W 


dU: 

aijn;  — ~  r 
^  ^  dx 


The  strain  energy  density  for  a  power-law  material  is  expressed  as 
(Kanninen  and  Popelar,  1985) 


W  = 


N 

N+1 


N+l 


(3.10) 


Equations  (3.9)  and  (3.10)  lead  to  the  following  general  representation  of 
stress,  strain  and  displacement  fields  near  the  crack  tip  (Rice  and 
Rosengren,  1968): 


aij(r,0;a,N)  =  OYK(a,N) 
a,(r,a;a,N)  =  aYK(a,N)  I,(0;N) 

eij(r,0;a.N)  =  aeYK^(a,N)  £..(0.^) 

Ui(r,0;a,N)  =  aeYK^(a,N)  Uij(0;N) 


(3.11) 


where  the  coordinate  functions,  Zij,L^,Eij,U j,  and  the  amplitude  of  the 
dominant  singularity,  K,  are  to  be  determined  by  solving  a  boundary  value 
problem.  Note  that  the  same  result  on  the  dominant  field  singularity  was 
obtained  by  Hutchinson  (1968a)  when  analysing  the  first  eigenvalue  of 
corresponding  differential  equation  for  the  Airy  function. 


3.2.  Formulation  of  the  Boundary  Value  Problem 


For  further  stress  and  strain  analysis,  the  Airy  stress  function  T  is 
introduced  as  follows  (Hutchinson,  1968a;  Rice  and  Rosengren,  1968) 


'P(r,0)  =  Kr*aY<l)(e)  (3.12) 

The  stress  components  which  satisfy  the  equilibrium  equations  are 

_  I  9'*'  .  1  3"^  .  -  9  Cl  3 

^ee  ~  -  2  '  ®rr  -  ^  +2  .,^2  •  ^re  “  " 


ar^  r  ar 

or,  using  eqs  (3.11)  and  (3.12), 


r"  ae^ 


ar 


r  ae 


Eee  =  s(s-l)(|)  ,  =  s<J)  +  r  .  =  (l'S)<t)' 


(3.13) 


with  s=(2N+l)/(N+l). 

Substituting  stress  components  (3.13)  into  the  compatibility  equation 
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one  can  find 


V 


d9^  N+). 


'e  '-'Try  -)  — e 

(N+l)" 


where  Sij  are  dimensionless  deviatoric  stress  components, 

Sij  =  ^ij  - 

Note  that  expressions  for  Eg  and  Sjj  are  different  for  plane  strain  and  plane 
stress.  After  certain  algebraic  rearrangements,  one  can  find  for  the  case  of 
plane  strain 


|<I„-Iesr+  31,2 


1/2 


(N+1)^ 


3N‘ 


(N+l)^ 


-(f)' 


1/2 


and 


(3.15) 


de- 


N(N+2) 


(N+1)^J 


'vN-1/^.k'-.  2N+1  ,  4N  d  /vN-u.N  _ 

|<tJ  + - 2"**  r - 

L  t  (N+1)"  JJ  (N+l)^deV  2 


=  0  (3.16) 


Calculations  for  the  case  of  plane  stress  yield 
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r."-,  NjaN+ir  ^2  _  N(^tl)  ^ 


(3.17) 
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(N+l)^ 


(N+1) 
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N+1 
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3N' 


.^2 
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rd2 


N  > 


d02  N+1 


lN-Y2^"+(N±2)<M±i)^ 
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N 


(N+1) 
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N-1 
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(N+1)^ 


(N+1)^  d0 


(3.18) 
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Equations  (3.16)  and  (3.18)  are  forth-order  differential  equations  for 

the  function  resolved  for  the  fourth-order  derivative,  d'*({)/d0'^. 

The  boundary  conditions  for  the  stress  function  can  be  formulated 
from  the  following  considerations.  Since  the  crack  surfaces  are  traction- 
free,  Oe0(+7c)='Cre(±7t)=O,  or  making  use  of  eqs  (3.13),  we  get  in  terms  of  the 
stress  function 

(j)(±tc)  =  (j>'(±rc)  =  O  (3.19) 

It  is  evident  that  the  stress  function  is  to  be  an  even  function  of  0  for  Mode 
I  crack  problem  and  an  odd  function  of  0  for  Mode  II  crack  problem.  Thus 
the  boundary  conditions  are 

<t)'(0)  =:  <j)"'(0)  =  4>(7C)  =  ^'(%)  =  0  (3.20) 

for  Mode  I,  and 

<})(0)  =  (t>’'(0)  =  <t)(ic)  =  i^\n)  =  0  (3.21) 

for  Mode  II. 


3.3.  Calculation  of  the  Amplitude  of  the  HRR-Fields 

Since  the  boundary  value  problem  formulated  in  the  previous  section  is 
solved,  the  amplitude  K  of  the  dominant  singular  fields  (3.11)  is  the  only 
unknown  parameter.  It  was  shown  by  Hutchinson  (1968a)  that  the  K- value 
is  expressed  via  the  y-integral,  the  stress  function  and  its  derivatives,  as 
well  as  parameters  of  the  Ramberg-Osgood  equation.  Substituting  eqs 
(3.10)  and  (3.11)  into  eq.(3.8)  one  can  find 


J  =  aavEyKNIK 


N+l 


or  K  = 


aavEyKN)] 


1/(N+1) 


with 


«N)=  I  {  +  cose 

-  IZ„(U9-U',)  -  I,8(U,+U’e)l  sine  }  M 


(3.22) 


(3.23) 


The  I-value  depends  on  the  strain  hardening  exponent  N,  as  well  as  on  the 
deformation  mode  and  plane  strain/plane  stress  conditions.  The 
displacement  functions  and  their  derivatives  obtained  from  Cauchy's 
equations  are 


U,  =  J(N+1)I^ 


N-l 


...  .  2N-I-1  . 
^  (N+l)2 


U',  =  ^(N+D-^d^-i) 

4  de  ' 


...  .  2N+1  ^ 


+  f  (N+1)S^ 
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N-l 


.  2N+1 

+ - r  <p 

(N+l)2  ^ 


(3.24) 


Ue  =  ^  U'r  +  3XN-1  (f-  , 
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for  plane  strain,  and 

Ur  =fN+i)i:N-> 


^  (N±2)(2^ 
2(N+l)2 


U'r  =(N+1)-^(Z^'-^) 

de 


+  (N+1) 


(t>'"  + 


2(N+1)2 
(N+2)(2N+1)  ,, 


2(N+1)2 


Ufi  =  ^  U’r  +  3X^-1  , 


N 

2 


U'fl  = 


LN-1)_(_2N±1)  _ 

(N+l)2 


-Ur 


for  plane  stress. 

3.4.  Numerical  Solution  of  the  Boundary  Value  Problem 


(3.25) 


As  it  follows  from  the  previous  analysis,  derivation  of  the  dominant 
singular  terms  of  elastic-plastic  crack  tip  fields  is  reduced  to  determination 
of  the  stress  function  which  is  to  satisfy  the  fourth-order  differential 
equation  (3.16)  or  (3.18)  and  boundary  conditions  (3.20)  or  (3.21) 
(depending  upon  the  mode  of  deformation  and  plane  strain/plane  stress 
cond'tions).  Thus,  the  two-point  boundary  value  problem  is  to  be  solved. 
For  this  purpose,  the  fourth-order  Runge-Kutta  integration  scheme  and  the 

shooting  method  are  employed.  The  latter  allows  for  reducing  the  two- 
point  boundary  value  problem  to  Cauchy's  problem. 

For  Mode  II  crack  the  boundary  conditions  (3.21)  can  be 

reformulated  into  the  initial  conditions  as  follows 

<|)(0)  =  0;  <{)'(0)  =  (pf.  <{>"(0)  =  0;  0"'(O)  =  (P2  (3-26) 

where  cpi  and  (P2  are  determined  through  an  iteration  process  to  satisfy  the 
boundary  conditions  at  0=7r; 

(t)(7t)  =  (j)'(7t)  =  0  (3.19a) 

Here  two  notions  are  to  be  remarked: 

(1)  Only  cpi  to  (p2  ratio  affects  the  solution  of  the  boundary  value 
problem.  Indeed,  it  follows  from  eqs  (3.7),(3.1 1),(3,13),(3.15),(3.17),(3.22)- 
(3.25)  that 

Zij  ~  (pi,  Sg  ~  (pi.  Eij  "  (p„  Ui  ~  (pf .  I(N)  ~  (pf-"‘ ,  K~  cpi^  , 

and  no  absolute  values  of  stress,  strain  and  displacement  components 

depend  on  (pj.  Thus,  an  arbitrary  value  of  <Pi(or  92)  choosen.  In 

the  present  analysis,  by  analogy  with  Hutchinson  (1968a,b)  and  Shih 
(1974),  the  (pi-value  is  put  to  satisfy  the  normalization  condition 
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while  the  (p  2 -value  is  found  by  the  shooting  method  to  satisfy  the 
boundary  conditions  (3.19a). 

(2)  Both  boundary  conditions  (3.19a)  are  satisfied  simultaneously. 
The  proof  of  this  fact  is  given  by  Rice  and  Rosengren  (1968)  and  is  based 
on  the  proof  of  the  7-integral  path-independence  (Rice,  1968a).  Since  the 
crack  surfaces  are  traction-free,  the  following  relation  is  valid  at  Q=k  and 
for  any  r-value; 

«  3U0  n 

-Oee  -r —  +  T^re— —  =  0  , 
dr  dr 

and  "the  vanishing  of  Tre  on  the  crack  surface  implies  the  vanishing  of  O90" 
(Rice  and  Rosengren,  1968). 

In  solving  the  boundary  value  problem,  no  convergence  was  found  for 
certain  initial  values  of  (p  2  with  the  derivative  (t>"‘  tending  to  infinity.  Thus, 
an  appropriate  first  approximation  function  (p2=(P2(N),  both  for  plane  strain 
and  plane  stress,  has  been  derived  in  numerical  experiment  and  employed 
after  in  the  iterative  process.  The  iterative  process  is  stopped  when  the 
following  condition  is  satisfied 

max  {|({)(;c)|  ;  l((>’(ic)l}  <  5  . 

with  5=10-'*.  A  step  size  of  0.45°  is  found  to  be  the  best  choise  in  the 
Runge-Kutta  integration  procedure.  Further  diminishing  of  the  step  size 
does  not  affect  accuracy  of  the  solution.  Under  above  conditions,  the 
convergence  of  numerical  procedure  is  obtained  in  two  to  nine  iterations 
within  the  range  of  the  strain  hardening  exponent  1<N<80,  both  for  plane 
strain  and  plane  stress. 

The  variation  of  stress  and  strain  fields,  5)ij  and  Eij,  at  a  crack  tip  in 
Mode  II  for  plane  strain  and  plane  stress  conditions  are  presented  in  Fig. 8 
and  Fig.9  respectively  (in  those  figures  the  index  “t”  corresponds  to  “9”); 
computations  were  performed  for  the  range  of  the  strain  hardening 
exponent  from  N=1.01  to  N=80,  and  the  limiting  cases  may  be  refered 
respectively  as  to  purely  elastic  and  perfectly  plastic  solutions. 
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4.  Rate  Sensitivity  of  Material  Properties 

Three  materials  with  different  rate  sensitivity  of  their  plasticity  properties 
have  been  analysed,  and  the  constants  in  the  Ramberg- Osgood  equation 
(TY,a,N)  change  accordingly  as  functions  of  strain  rate: 

(1)  mild  steel',  the  basic  data  on  the  rate  sensitivity  of  the  yield 
stress  are  available  from  the  work  of  Campbell  and  Ferguson  (1970)  where 
steel  of  the  following  chemical  composition  was  studied:  0.12%C,  0.10%Si, 
0.62%Mn.  0.029%S,  <0.{)04%Pb.  Also  data  of  Klepaczko  (1969)  on  technically 
pure  iron  (0.05%C)  and  Tanaka  and  Kinoshita  (1967)  (0.03%C)  were  used  to 
study  strain  rate  effect  on  parameters  describing  stress-strain  relation; 

(2)  /«//>■  annealed  titanium  alloy  Ti6A14V;  the  data  on  behavior  of 
stress-strain  curves  at  different  strain  rates  are  available  from  the  paper 
by  Maiden  and  Green  (1966); 

(3)  aluminum  alloy  A1  6061-T6  (l%Mg,  0.6%Si,  0.25%Cu,  0.25%Cr); 
the  strain  rate  effect  on  the  yield  stress  and  stress-strain  curves  for  this 
alloy  was  studied  by  Maiden  and  Green  (1966),  some  additional  results  can 
be  found  in  the  work  of  Jiang  and  Chen  (1974). 

The  behavior  of  those  materials  has  been  analysed  for  the  case  of 

“4  3-1 

room  temperature  and  strain  rates  from  10  to  10  s  . 

Figure  1 1  shows  stress-strain  curves  for  technically  pure  iron 

obtained  in  pure  shear  at  the  strain  rate  range  7=  1.31x10-4. ..55  s- 1 
(Klepaczko,  1969).  For  this  material,  the  lower  (Tyl)  ^d  upper  (tyo)  yield 
stresses,  as  well  as  the  slope  of  stress-strain  curve,  are  changed 
significantly  with  increasing  strain  rate.  The  experimental  data  of 

Klepaczko  (1969)  on  variation  of  Tyl  versus  log  y  are  shown  in  Fig.r2a, 
corresponding  results  of  Campbell  and  Ferguson  (1970)  and  Tanaka  and 
Kinoshita  (1967)  are  also  presented.  Note  that  original  results  of  the  latter 
authors  obtained  in  compression  are  transformed  here  to  the  stress-strain 
diagram  in  shear  according  to 

Tyl  =  Oyl/'v/s'.  Y  =  \'3e  (4.1) 

An  averaged  Tyl  versus  log  y  curve  is  fitted  to  the  experimental  data  in 
Fig.  12a  within  the  strain  rate  range  from  10' 4  to  lO^  s'F 

Based  on  results  of  Klepaczko  (1969),  Fig. 11,  the  Ramberg-Osgood 
constants  N  and  a  in  eq.(3.1),  have  been  found  and  they  are  shown  versus 
strain  rate  in  Fig.  12b.  It  was  observed  in  many  studies  on  mild  steel  (see, 
e.g.,  Klepaczko,  1969)  that  with  increasing  strain  rate  the  slope  of  the  strain 
hardening  part  of  the  stress-strain  curve  diminishes.  Above  certain  level  of 

strain  rate  ('jp’lO^s*^)  when  adiabatic  heating  is  of  importance,  this  part  of 
the  stress-strain  diagram  shows  a  negative  slope.  Since  it  is  impossible  in 
the  present  analysis  to  model  "negative"  strain  hardening,  it  is  assumed 

that  the  value  of  N  tends  to  infinity  at  high  strain  rates  (V>102  5-I).  Xhis 
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case  corresponds  to  the  model  of  elastic-perfectly*plastic  material.  With 

this  assumption  taken  into  account,  the  N  versus  log  y  curve  in  Fig.  12b  is 

extrapolated  to  y=10^  S'^.  An  appropriate  extrapolation  of  the  constant  a 

versus  log  y  is  also  shown.  Based  on  results  presented  in  Fig.l2a,b,  one  can 
conclude  that  there  is  strong  effect  of  strain  rate  on  the  Ramberg-Osgood 
parameters  for  mild  steel.  Both  the  lower  yield  stress,  Xy,  and  strain 

hardening  exponent,  N,  are  rapidly  increasing  with  log  y,  while  the  value  of 
a  is  slightly  decreasing. 

Figure  13  shows  test  data  of  Maiden  and  Green  (1966)  obtained  in 

compression  for  titanium  alloy  Ti6A14V  in  the  range  of  strain  rate  e  = 

4x10*3. .,20  s*l,  or  after  transformation  according  to  (4.1),  y-  6.9x10*3. ..34. 6 
s'F  The  strain  rate  dependence  of  the  yield  stress,  Xy,  for  this  alloy  is 

shown  in  Fig. 14a.  Here  the  yield  stress  in  shear  is  found  as  Xy  =  CT0.2W  3  . 

The  Xy  versus  log  y  curve  is  extrapolated  to  two  extreme  limits  Y=10*^s*i 

and  Y=103  s-F  Variation  of  the  constants  in  the  Ramberg-Osgood  equation 
is  presented  in  Fig.  14b.  One  can  conclude  that  no  much  effect  of  strain  rate 

on  a  and  N  is  observed.  Both  the  N  versus  log  y  and  a  versus  log  y  curves 
are  approximated  by  straight  lines  and,  thus,  extrapolated  within  the 
whole  range  of  the  strain  rate  considered  (Fig.  14b). 

Test  data  on  Al  6061-T6  (Maiden  and  Green,  1966)  in  compression  at 

the  range  of  strain  rate  £=9x10*3. ..910  s**  or,  equivalently,  Y=l. 56x10*2... 
1.58x103  s*l  (Fig. 15),  reveal  no  strain  rate  effect  on  the  stress-strain  curve. 
The  curve  is  approximated  by  the  Ramberg-Osgood  equation  (3.1)  with 
Xy  =  175  MPa,  N=4.3  and  a  =  1.54.  Similar  results  but  with  a  small  rate 
sensitivity  of  the  yield  stress,  O  y ,  have  been  reported  for  the  same 

material  by  Jiang  and  Chen  (1974),  within  the  range  of  strain  rate 

e=l .6x10*3. ..210  s^  or y  =2. 8x10*3. ..360  s*F  Their  data  are  presented  in 

Fig.  16  as  Xy  versus  log  y  and  show  a  linear  dependence  of  the  yield  stress 
versus  logarithm  of  strain  rate. 

It  is  important  to  note  that  in  Fig.  15  the  engineering  stress-strain 
curves  which  are  presented  with  no  account  for  adiabatic  heating  at  high 
strain  rates.  In  fact,  if  the  effect  of  adiabatic  heating  is  considered,  an 
aluminum  alloy  reveals  strain  rate  sensitivity  of  its  mechanical  properties. 

Note  also  that  at  high  strain  rates  (  y  >103  s**),  the  0.2%  proof  stress  of 
aluminum  alloy  becomes  rate  sensitive  and  its  behavior  is  similar  to  that  of 
yield  stress  observed  for  steels. 
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5.  Strain  Rate  Effect  on  Evolution  of  the  Crack  Tip 
Plastic  Zone  for  Structural  Materials 
Via  SSY  Approximation 

In  the  previous  chapters  rate  sensitivity  of  mechanical  properties  of  mild 

steel,  titanium  and  aluminum  alloys  in  the  strain  rate  range  7=  10"*...  103  yi 
has  been  discussed,  as  well  as  the  two-parameter  model  of  the  plastic  zone 
near  a  crack  tip  in  Mode  II  has  been  analyzed  within  the  SSY 
approximation.  In  that  model  both  size  and  shape  of  the  plastic  zone 
depend  strongly  on  the  ratio  of  the  applied  shear  stress  to  the  yield  stress, 
'Cappi^'^Y'  as  well  as  on  the  crack  geometry  specified  by  the  value  of  the 
stress  intensity  factor,  Ku.  In  the  next  part  the  analysis  of  evolution  of  the 
plastic  zone  near  the  crack  tip  in  Mode  II  with  parametric  variation  of  the 
strain  rate  is  performed  for  different  materials. 

The  crack  geometries,  which  are  of  practical  interest,  are:  (1)  semi¬ 
infinite  crack  in  an  infinite  plate  and  (2)  finite  length  edge  crack  in  a  strip. 
The  first  geometry  represents  the  limiting  case  of  a  very  large  crack  which 
is  typical  for  problems  arising  in  geomechanics,  pipeline  safety  assessment 
etc.  The  second  one  corresponds  to  specimen  geometries  used  in  fracture 
toughness  tests  in  Mode  II  (Banks-Sills  and  Arcan.  1986;  Banks-Sills  and 
Sherman,  1991;  Davies  et  ai,  1985;  Kalthoff  and  Winkler,  1987;  Maccagno 
and  Knott,  1992;  Mason  et  ai,  1992;  Nishioka  et  ai.  1991;  Tohgo  et  ai, 
1989;  Tohgo  and  Ishii,  1992).  The  geometry  used  by  Kalthoff  and  Winkler 
(1987)  and  Nishioka  et  al.  (1991)  with  crack  length  a  =  50  mm  and 
specimen  width  W  =  100  mm  is  studied  in  the  present  work. 

Analysis  of  the  evolution  of  plastic  zone  with  strain  rate  is  carried  out 

for  constant  value  of  applied  stress,  i.e.  Tappi(7)=const,  while  the  yield  stress 
as  a  function  of  the  strain  rate  is  involved.  Thus,  the  load  level  is 

characterized  by  the  ratio  of  tjppi  to  the  yield  stress  Tys  in  test  with  7=1  O'"* 
S‘f  which  is  refered  to  as  quasi-static  case.  Two  different  stress  levels  are 
analysed:  Tappi/XYs=0.2  and  0.5. 

To  characterize  the  evolution  of  plastic  zone,  three  parameters  are 
introduced: 

(1)  plastic  zone  size  ahead  of  the  crack  tip,  rp(0=O),  which  is  normally 
the  maximum  size; 

(2)  half  the  plastic  zone  hight,  hmax  (since  the  plastic  zone  is  symmetric 
with  respect  to  the  crack  line,  its  total  maximum  hight  equals  to  2hinax); 

(3)  the  ratio  of  the  plastic  zone  size  behind  the  crack  tip  to  that  ahead 
of  the  crack  tip,  rp(0=7c)/rp(9=O);  this  parameter  characterizes  asymmetry  of 
the  shape  of  plastic  zone. 

In  the  case  of  a  semi-infinite  crack  the  values  of  rp(0)  and  h^ax  are 
normalized  by  the  crack  length,  a.  Analyses  are  carried  out  for  plane  strain 
conditions  with  Poisson's  ratio  7=0.33. 
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Estimations  of  the  strain  rate  effect  on  the  evolution  of  plastic  zone 
near  the  tip  of  a  semi-infinite  and  an  edge  crack  in  mild  steel  are 
presented  in  Fig.  17  and  Fig.  1 8.  One  can  find  significant  reduction  of  the  size 
of  plastic  zone  with  increasing  strain  rate  from  lO"'*  to  10^  s'l.  For  the 
stress  level  T3ppi/tYS=0.2  (Figs  17a  and  18a),  revluction  of  rp(0)  and  hmax. 

within  the  range  of  Y  considered,  reaches  value  of  five  times  for  the  two 
crack  geometries  studied,  while  the  ratio  rp(ji)/rp(0)  is  kept  nearly  constant 
and  equal  1,  i.e.  the  plastic  zone  is  nearly  symmetric  with  respect  to  the  y- 
axis  (0=3i/2).  For  the  stress  level  iappi/tYs=0-5  (Figs  17b  and  18b),  there  is 
stronger  rate  effect  both  on  the  plastic  zone  size  and  shape.  If  the  strain 
rate  is  increased  from  10-'*  to  10^  s-^  the  values  of  rp(0)  and  hmax  for  a 
semi-infinite  crack  are  diminished  respectively  7.4  and  6.0  times.  For  an 
edge  crack  this  reduction  reaches  5.6  and  5.2  times.  One  can  see  also  from 
Figs  17b  and  18b  a  significant  change  in  the  shape  of  plastic  zone  which  is 
asymmetrical  in  quasi-static  loading  and  nearly  symmetrical  at  high  strain 
rate.  The  asymmetry  of  the  plastic  zone  for  quasi-static  loading  is 
characterized  by  factor  rp(7c)/rp(0)=0.47  for  a  semi-infinite  crack  and  0.76 
for  an  edge  crack. 

Significant  strain  rate  effect  on  the  size  and  shape  of  plastic  zone 
found  for  mild  steel  is  a  consequence  of  strong  rate  sensitivity  of  the  lower 

yield  stress.  The  latter  is  increased  2.2  times  with  increasing  7  from  10-*to 
10^  S'*.  For  titanium  alloy  corresponding  increase  of  the  yield  stress 
reaches  1.4  which  results  in  a  lower  reduction  of  the  plastic  zone  in 
comparison  to  steel.  At  tapp!/XY=0.2,  the  size  of  plastic  zone  is  diminished 
twice  for  both  crack  configurations  (Figs  19a  and  20a).  At  Tappi/XY=0.5,  the 
value  of  rp(0)  is  decreased  respectively  2.6  and  2.2  times  for  a  semi- 
infinite  and  an  edge  crack,  while  the  value  of  hmax  is  diminished  twice  for 
both  crack  configurations  (Figs  19b  and  20b).  The  change  in  the  factor 
rp(7c)/rp(0)  is  similar  to  that  observed  for  mild  steel  with  the  shape  of 
plastic  zone  being  more  asymmetric  at  higher  strain  rates  and  Tappi/XY=0.5. 

As  it  is  mentioned  in  Chapter  4,  no  rate  sensitivity  of  the  yield  stress 

for  aluminum  alloy  is  found  at  y  =10-**...  10^  s-i.  Thus,  one  can  expect  no 
strain  rate  effect  on  the  size  and  shape  of  plastic  zone  within  this  range  of 

y.  The  predicted  values  of  the  parameters  of  the  plastic  zone  geometry  are 
listed  below; 
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'‘app!''ry— 0.2 

tappi/'tY— O.O 

Semi-infinite  crack 

rp(0)/a 

0.017 

0.168 

htnax/ <2 

0.008 

0,045 

rp(a)/rp(0) 

0.929 

0.470 

Edge  crack 

rp(0),  mm 

1.605 

11.39 

bniax»  fnm 

0.730 

4.64 

rp(rc)/rp(0) 

0.994 

0.79 

Note  that  the  present  approach  to  estimate  the  size  and  shape  of 
plastic  zone  takes  no  account  of  stress  redistribution  due  to  plasticity.  If 
Irwin’s  correction  is  applied,  the  value  of  rp  ahead  of  the  crack  tip  is  twice 
larger  then  predicted  by  Figs  17  to  20.  while  the  rp(:t )  and  hmax  values 
may  remain  unchanged.  Thus,  more  significant  plastic  zone  asymmetry  can 
be  expected  with  the  plastic  zone  being  predominantly  concentrated  ahead 
of  the  crack  tip.  For  instance,  with  Irwin’s  correction,  the  rpfir)/rp(0)  factor 
for  an  edge  crack  in  mild  steel  at  tappi/Tys=0.5  reaches  value  0.38.  This 
result  is  in  qualitative  agreement  with  the  finite-element  solution  of 
Banks-Sills  and  Sherman  (1990).  (Since  a  specimen  w'ith  stress 
concentrator  is  analyzed  by  those  authors,  quantitative  comparison  of 
results  is  not  available.) 

It  seems  reasonable  to  say,  by  analogy  with  Mode  I  crack,  eqs  (1.4) 
and  (1.5),  that  cenain  relations  between  the  maximum  size  of  plastic  zone 
(>"p,max)-  material  fracture  toughness  (/fiic)  and  minimum  specimen 
thickness  (B)  must  exist.  If  it  is  the  case  and  the  ratio  rp(0)/B  is  considered 
as  a  measure  of  plasticity  effects  at  the  crack  tip,  the  results  on  variation  of 
rp  with  strain  rate  can  be  used  in  predicting  the  minimum  thickness  (or 

width)  of  specimen  versus  7.  For  example,  at  the  stress  level  Xappi/'CYS=0.5 
(which  is  normal  in  Mode  II  fracture  toughness  tests;  see  Banks-Sills  and 
Sherman.  1991),  our  results  for  an  edge  crack  in  mild  steel  (Fig.  18b) 

predict  the  specimen  thickness  at  high  strain  rate  (7=10^S'^)  5.6  times 

smaller  than  in  quasi-static  test  (7=10 For  titanium  alloy  at  the  same 
conditions  (Fig.20b),  reduction  of  the  specimen  size  is  2.2  times. 
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6.  Strain  Rate  Effect  on  Behavior  of  the  Crack  Tip 
Fields  in  Mode  II  Via  the  HRR-Solution 

T  ie  results  on  evolution  of  the  plastic  zone  with  strain  rate  as  a  parameter 
pies  anted  in  Chapter  5  are  based  exclusively  on  the  rate  sensitivity  of  the 
yiel<i  stress.  Some  additional  information  on  the  behavior  of  crack  tip  fields 
at  different  strain  rates  can  be  obtained  from  analysis  of  the  HRR-solution, 
eqs  (3.11),  which  incorporate  parameters  of  the  stress-strain  curve. 
However,  it  is  important  to  note  that  the  HRR-solution  cannot  be  directly 
employed  lo  evaluate  the  size  of  the  plastic  zone  near  crack  tip.  Indeed,  the 
plastic  zone  boundary  is  interpreted  as  a  line  of  the  constant  equivalent 
stress,  i.e.  with  a|*=0.  Since  the  elastic  stress  and  strain  components 

have  been  neglected  when  deriving  the  HRR-field  (see  eq.(3.7)),  the  latter 
seems  to  be  valid  only  at  high  stress  and  strain  levels, 

cjg  >  Oy  and  eP  »  Cy  , 

or,  equivalently,  in  a  small  vicinity  of  the  crack  tip  with  a  radius  bein^, 
much  smaller  the  plastic  zone  size,  i.e.  rjjRj^  «  rp.  Moreover,  eqs  (3.11)  yield 
only  the  dominant  components  of  the  crack  tip  field.  The  range  of 
dominance  of  the  KRR-solution  is  often  determined  as  (e.g.,  Du  et  al.,  1991) 


with  P  ranging  from  about  2  to  5.  In  the  case  of  numerical  example 
presented  in  Fig.  18b  (Mode  II  crack  of  finite  length  in  mild  steel)  one  can 
find  from  eq.(6.1)  Thrr  =  0.05. ..0.1  mm. 

Rice  and  Rosengren  (1968)  calculated  boundaries  r(,(0)  of  the  equal 
equivalent  stress  versus  the  strain  hardening  exponent  in  Mode  I.  Their 
results  on  the  rg-evolution  are  in  qualitative  agreement  with  direct 
computations  of  the  size  and  shape  of  plastic  zone  near  the  crack  tip  (Shih, 
1974).  Thus,  a  qualitative  analysis  of  the  behavior  of  plastic  zone  near  the 
crack  tip  at  different  strain  rates  can  be  performed  employing  the  HRR- 
solution.  However,  it  seems  to  be  more  rational  in  such  analysis  to 
characterize  me  evolution  of  crack  tip  fields  via  the  equivalent  plastic 
strain  which  simplifies  identification  of  a  corresponding  point  of  stress- 
strain  curves  obtained  at  different  strain  rates. 

6.1.  Evolution  of  the  Strain  Field  near  a  Crack  Tip  in  Mode  II 
for  Mild  Steel  and  Titanium  Alloy 

Numerical  examples  of  the  strain  field  near  a  crack  tip  are  given  for  the 
case  of  a  finite  length  edge  crack  in  a  strip  with  a  =  50  mm  and  W  =  100 
mm.  Corresponding  assessments  of  evolution  of  the  plastic  zone  founded  on 
the  SSY-approximation  are  presented  in  Figs  18  and  20  for  mild  steel  and 
titanium  alloy  respectively. 
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Making  use  of  eqs  (3.6),(3.1 1).(3-22)  and  (2.7),  one  can  find  the 
radius  of  a  boundary  of  the  equal  equivalent  plastic  strain  ej  =  (with 
P£=const)  as  follows 


1  -  v' 


3aI(N) 


a  \n+i)/n  ^n+1 


a/3  ^ 


'^n> 


(6.2) 


Here  the  small  scale  yielding  and  plane  strain  conditions  at  the  crack  tip 
are  assumed;  also  the  Huber-Mises  yield  criterion  is  employed. 

Figure  21  presents  variation  of  rg-boundaries  at  the  crack  tip  in  mild 
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steel  at  the  strain  rate  range  1.31x10  ...10  s’  .  Two  deformation  levels 
are  ane’ysed:  =  0.03  and  0.06  which  correspond  approximately  to  points 

with  'f  =  0.05  and  0.1  on  the  stress-strain  diagrams  in  Fig  ll.  The  results  in 
Fig. 21  demonstrate  more  complicated  and  somewhat  different  behavior  of 
the  crack  tip  fields  versus  the  strain  rate  than  those  presented  in  Fig.l8a,b. 
With  increasing  strain  rate  the  dimensionless  radius 


of  the  tg-boundaries  ahead  of  the  crack  tip  sligtly  diminishes  and  after  takes 


.0  .\ 

the  maximum  value  at  y=  2.85x10  "  s  .  The  values  of  Rg,ina)c  0.0072  and 
0.003  for  £p  =  0.03  and  0.06  respectively.  Further  increase  in  the  strain  rate 
results  in  fast  diminishing  of  the  Rg  value  ahead  of  the  crack  tip.  This 
nonmonotonic  variation  of  the  Rg^^iax  value  is  a  consequance  of  a  complex 
dependence  of  plasticity  parameters  (Oy  ,  a,  N)  upon  the  strain  rate  (Fig.  12). 

At  lower  strain  rates  (1.31x10”*  s’*  to  2.85x10'^  s  ^)  the  value  of  a  is 
nearly  constant,  and  increase  of  Oy  and  N  with  strain  rate  leads  to  Rg.max  6e 
almost  unchanged.  At  higher  strain  rates  the  HRR-solution  predicts 
substantial  reduction  in  the  value  of  Rg.max  due  to  decreasing  of  a  and 
increasing  of  Oy,  which  effect  is  much  stronger  than  expansion  of  the  Rg- 
zone  due  to  increase  of  the  strain  hardening  exponent. 

In  contrast,  the  Rg  value  behind  the  crack  tip  and  the  maximum  hight 


of  tg-boundaries  decrease  monotonically  with  increasing  y  .  Another 
important  feature  of  tg-boundaries  behavior  consists  in  fast  diminishing  of 

the  rg(:c)/rg(0)  ratio  which  changes  from  0.24  at  y=  1.31x10'“*  s'^  to  0.05  at 

7=  10  s  .  Thus,  if  a  certain  correlation  between  and  tp  values  exists,  one 
can  conclude  that  at  high  strain  rates  almost  the  whole  plastic  zone  area  is 
concentrated  ahead  of  the  crack  tip;  note  that  the  SSY-approximation 
(Fig.  18)  predicts  the  plastic  zone  at  high  strain  rates  nearly  symmetrical 
with  respect  to  the  y-axis. 
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Numerical  results  on  evolution  of  Tj-boundaries  at  the  crack  tip  in 
titanium  alloy  Ti6A14V  are  presented  in  Fig.22  for  e?  =  0.02.  This 

deformation  level  corresponds  to  =  0.035  and  total  deformation  - 
0.032. ..0.036  in  Fig.  13.  Comparing  to  mild  steel,  no  significant  rate 
sensitivity  of  r^-boundaries  at  the  crack  tip  in  titanium  alloy  is  observed: 

with  increasing  strain  rate  from  10  s  to  10  s  ,  the  maximum  size  and 
hight  of  the  tg-boundary  is  diminished  approximately  15%,  while  the  size  of 
the  tg-boundary  behind  the  crack  tip  is  decreased  twice.  This  is  due  to 
relatively  weak  rate  sensitivity  of  the  plasticity  parameters  for  titanium 
alloy. 

6.2.  Comparison  of  Approximate  Results  on  the  Plastic  Zone 
Analysis  with  Numerical  Solutions 


Only  a  few  numerical  solutions  on  the  plastic  zone  near  a  crack  tip  in  Mode 
II  are  known.  Shih  (1974)  and  Pan  and  Shih  (1990)  published  results  on 
variation  of  the  plastic  zone  geometry  versus  the  strain  hardening  exponent 
at  constant  a  and  Oy.  Their  finite-element  solutions  are  constructed  within 
the  framework  of  the  boundary  iayer  formulation  for  small  scale  yielding 
and  plane  strain.  Banks-Sills  and  Sherman  (1990)  obtained  an  elastic-plastic 
stress  field  for  the  specimen  shown  in  Fig. 6. 17.  Different  levels  cf  the 
applied  load  were  analysed;  Xg^ppj/Xy  =  0.35. ..0.52,  where  x^ppj  is  a  nominal 
stress  defined  as 


t  is  the  specimen  thickness,  W  is  the  specimen  width  in  a  cracked  section 
(W  =  30  nun  in  Fig.6.17).  The  lower  and  upper  load  levels  may  be  refered 
respectively  to  as  a  small  scale  yielding  and  fully  plastic  solutions. 

Plastic  zone  boundaries  versus  N  are  shown  in  Fig. 23,  after  Shih 
(1974)  and  Pan  and  Shih  (1990).  In  both  papers  the  constant  a  in  the 
Ramberg-Osgood  equation  is  equal  3/7.  Variation  of  the  plastic  zone  versus 
the  load  level  is  presented  in  Fig.24,  after  Banks-Sills  and  Sherman  (1990). 
Their  results  are  obtained  for  aluminum  alloy  A1  7075-T7351  with 
parameters  of  the  Ramberg-Osgood  equation  found  to  be  a=1.25,  N=14. 
Table  6.1  lists  dimensionless  parameters  of  the  plastic  zone  geometry 


^p*max 


=  r„(0) 


iCn 


Hp=hp 


and  ro(7t)/r„(0) 


found  from  the  papers  mentioned  above,  as  well  as  results  obtained  in 
Chapter  2  founded  on  the  two-parameter  crack  tip  stress  field 
representation  (Fig.5b)  and  with  Irwin's  plasticity  correction.  One  can 
conclude  from  analysis  of  numerical  solutions  that  both  effects  of  the  strain 
hardening  exponent  and  load  level  on  the  plastic  zone  behavior  are  similar: 
increase  of  N  or  t^ppi/xy  yields  substantial  growth  of  the  plastic  zone  size 


Parameters  of  the  plastic  zone  geometry 


Table  6,1 


Reference  t 

appl'^’^Y 

a 

N 

Hp 

rp(7t)/rp(0) 

Shih  (1974) 

3/7 

1 

0.16 

0.067 

1.0 

- 

3 

0.24 

0.080 

0.5 

- 

13 

0.29 

0.076 

0.21 

Pan  and  Shih 

3/7 

3 

0.17 

0.064 

0.76 

(1990) 

10 

0.21 

0.061 

0.32 

Banks-Sills  and 

0.35 

1.25 

14 

0.48 

0.11 

0.24 

Sherman  (1990) 

0.45 

0.69 

0,12 

0.15 

0.48 

0.54 

0.13 

0.13 

0.52 

0.43 

0,10 

0.15 

Present  study, 

0.2 

0.34 

0.071 

0.46 

2 -parametrical 

0.3 

0.37 

0.074 

0.39 

elastic  stress  field 

0.5 

0.53 

0.085 

0.23 

representation 

Present  study, 

3/7 

** 

\ 

0.26 

0,075 

0.28 

tg-boundaries 

- 

13 

0.30 

0.067 

0.07 

according  to 

eq.(6.3) 

Present  study, 

2n 

3 

0.24 

0.072 

0.29 

rj-boundaries 

13 

0.38 

0.091 

0.06 

ate|’/eY=  1 

ahead  of  the  crack 

tip  and 

diminishing 

of  the  rp( 

:t)/rp(0) 

ratio.  At  the 

time  only  slight  change  in  the  Hp  is  observed. 

One  can  find  certain  disagreement  in  reviewed  numerical  results.  The 
data  of  Shih  (1974)  and  Pan  and  Shih  (1990)  for  identical  conditions  (a=3/7, 
N=3)  differ  up  to  50%  when  analysis  is  made  for  the  Rp,max  value  and 
rp(rt)/rp(0)  ratio.  Comparison  of  the  plastic  zone  solution  of  Banks-Sills  and 
Sherman  (1990)  with  results  cf  Shih  (1974)  and  Pan  and  Shih  (1990)  is  not 
available  because  of  the  difference  in  values  of  a  used  in  those  analyses. 
Nevertheless,  the  asymmetry  factor,  rp(7c)/rp(0),  in  the  SSY  analysis 
('Cjppj/Xy=0.35)  of  Banks-Sills  and  Sherman  (1990)  is  very  close  to  that  in 
Shih’s  (1974)  solution  with  N=13. 

It  should  be  noted  befor  discussing  results  of  the  present  study  on  the 
plastic  zone  behavior  versus  the  load  level,  that  the  approach  employed  in 
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Chapter  2  together  with  Irwin's  plasticity  correction  is  equivalent  to  the  use 
of  the  elastic-perfectly-plastic  model  (N=<»).  This  is  close  to  the  case 
considered  by  Banks-Sills  and  Sherman  (1990),  and  comparison  of  our 
results  with  the  finite-element  solution  of  Banks-Sills  and  Sherman  (1990) 
seems  to  be  resonable.  One  can  find  from  the  data  presented  in  Table  6.1 
that  values  of  Rp^max  determined  in  the  present  study  are  within  20%  of  the 
numerical  results,  while  the  error  in  estimation  of  Hp  and  rp(7t)/rp(0)  is 
higher  (20. ..50%).  Thus,  for  materials  with  high  values  of  the  strain 
hardening  exponent  the  two-parameter  crack  tip  field  characterization 
together  with  Irwin's  plasticity  correction  yield  realistic  assessments  of  the 
plastic  zone  size  ahead  of  the  crack  tip  in  Mode  II. 

6.3.  Correlation  of  r^  and  r Boundaries  with  the  Plastic  Zone 
Geometry 

It  is  interesting  to  compare  boundaries  of  equal  equivalent  stress  and 
plastic  strain  with  plastic  zone  configurations.  Figure  25  shows  variation  of 
rQ-boundaries  versus  the  level  of  .  Calculations  were  carried  out  for 
a  =3/7  and  two  values  of  the  strain  hardening  exponent:  N=3  and  13.  After 
comparing  of  those  results  with  data  of  Shih  (1974)  (see  Fig.23a  and  Table 
6.1),  one  can  conclude  coincidence  of  tj,  and  tp  values  ahead  of  the  crack  tip 

at  Og  =  l.SOv  for  N=3  and  at  Og  =  l.OSoy  for  N=13.  It  is  evident  that  at  N->«> 
the  condition  r^  =  tp  holds  at  =  Oy*  fhus  for  N>3  the  following  relationship 

=  Oy  (I  +  0-9/N)  (6.3) 

seems  to  be  valid  to  estimate  the  rp(0)  value  from  the  HRR-field.  The 
numerical  results  on  configuration  of  tf^-boundaries  for  the  stress  level 
determined  by  eq.(6.3)  are  presented  in  Table  6.1.  Note  that  not  only  values 
of  ro(0)  but  also  the  maximum  hight  of  the  r„-boundaries  are  in  good 
agreement  with  ip-solution  of  Shih  (1974).  At  the  same  time  the  value  of 
r<j(ji)  is  underestimated.  This  is  a  consequance  of  the  fact  that  the  HRR- 
solution  involves  only  dominant  singular  terms  of  the  stress  field.  One  could 
expect  that  the  use  of  nonsingular  terms  in  eqs  (3.11)  would  lead  to  a  more 
realistic  assessments  of  as  well  as  tp-boundaries. 

For  the  same  conditions  (a=3/7,  N=3  and  13),  Fig.26  shows  variation  of 
fj-boundaries  versus  the  level  of  ej  .  One  can  see  a  good  agreement  between 
tg-values  ahead  of  the  crack  tip  and  tp-solution  of  Shih  (1974)  at  eP  =  (1  ... 
1.5)  Ey  (Table  6.1).  Since  the  condition  r^iO)  =  rp(0)  is  satisfied  at  nearly 
constant  level  of  eP  ,  the  use  of  tj-boundaries  to  characterize  the  plastic  zone 
behavior  versus  N  seems  to  be  preferable.  However,  this  conclusion  is  valid 
only  for  materials  which  stress-strain  curves  can  be  approximated  by  the 
Ramberg-Osgood  equation  at  low  deformation  levels,  eP  =  fiy  .  It  is  not  the 
case  for  mild  steel  which  stress-strain  diagram  in  low  strain  level  is 
nonmonotonic  (Fig.  11). 
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7,  Conclusions 

In  the  present  study  an  approximate  parametric  analysis  of  the  strain  rate 
effect  on  the  crack  tip  fields  in  Mode  11  is  carried  out.  Evolution  of  the 
plastic  zone  and  crack  tip  stress  and  strain  fields  at  different  strain  rates  is 
modelled  via  corresponding  quasi-static  solutions  which  couple  strain  rate 
sensitivity  with  plasticity  parameters  (Xy.  cc,  N).  An  attempt  is  made  to 
predict  the  size  and  shape  of  the  plastic  zone  within  the  SSY  approximation 
based  on  the  two-parameter  characterization  of  the  crack  tip  field,  eq. 
(2.16).  For  materials  with  a  high  strain  hardening,  such  approach  together 
with  Irwin's  plasticity  correction  is  found  to  provide  an  adequate 
evaluation  of  the  plastic  zone  size  ahead  of  the  crack  tip  as  well  as  the 
maximum  hight  of  plastic  zone.  Our  results  on  rp  j^ax'^ariation  versus  the 
load  level  are  within  20%  of  the  numerical  solution  by  Banks-Sills  and 
Sherman  (1990),  for  t^pp^/xy  ranging  from  0.35  to  0.52.  By  analogy  v/ith  the 
case  of  Mode  I  (Larsson  and  Carlsson,  1973),  the  two-parameter 
characterization  of  the  crack  tip  field  can  be  employed  in  finite-element 
analysis  through  the  boundary  layer  formulation.  In  this  way  more 
realistic  plastic  zone  assessments  can  be  obtained  with  minimum 
computational  efforts. 

Making  use  of  the  SSY  approximation,  evolution  of  the  plastic  zone 
near  the  crack  tip  is  analysed  at  different  strain  rate  for  materials  with 
different  rate  sensitivity  of  the  yield  stress.  Significant  decrease  of  the 
plastic  zone  size  with  increasing  strain  rate  is  found  for  mild  steel.  At  the 
load  level  Xappi/Xy  =  0.5,  which  is  typical  for  fracture  toughness  tests, 

I  3*1 

increase  of  the  strain  rate  from  10  s’  to  10  s*  leads  to  5,6...  7.4  times 
reduction  in  the  maximum  plastic  zone  size  (depending  upon  a  crack 
geometry).  The  rate  sensitivity  of  the  plastic  zone  size  is  stronger  at  higher 
load  levels.  For  titanium  alloy,  the  reduction  factor  of  the  plastic  zone  size 
is  2.2  at  the  same  range  of  the  strain  rate.  No  rate  sensitivity  of  the  plastic 
zone  is  predicted  for  the  aluminum  alloy.  Since  the  maximum  size  of  the 
plastic  zone  can  be  considered  as  a  measure  of  plasticity  effects  at  the 
crack  tip,  one  can  expect  reduction  of  the  minimum  specimen  thickness  for 
ifjl(n-measurement  at  increased  strain  rates  as  proportional  to  the  change  of 

^p.max- 

When  the  HRR-solution  is  employed  to  characterize  the  crack  tip 
fields,  evolution  of  boundaries  of  equal  equivalent  stress  and  plastic  strain 
can  be  derived.  For  mild  steel,  the  maximum  size  of  the  tg- boundary 
changes  in  nonmonotonic  way  versus  strain  rate  (Fig.21),  with  the 

maximum  value  at  y=  2.85x10'^  s'^  for  =  0.03  and  at  y  =  0.44  s’^  for  e|  = 

-2  -3 

0.06.  This  maximum  value  found  from  Fig.21  is  of  the  order  (10  ...  10  )a  , 
which  is  close  to  the  process  zone  size.  For  deformation  levels  considered 
decrease  of  the  tj-boundary  for  mild  steel  within  the  range  of  strain  rate 
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from  10  "^  s'^  to  10^  s'^  constitutes  40...  50%.  In  contrast,  the  size  of 
boundary  for  the  titanium  alloy  is  changed  monotonically  with  strain  rate, 
and  the  total  reduction  of  the  r^  max^value  within  the  range  of  strain  rate 

from  lO'"^  s*^  to  10^  s'^  is  10%. 

From  the  point  of  view  of  experimental  investigation  of  crack 
initiation  and  fracture  toughness  measurement  in  Mode  II,  application  of 
three  types  of  specimens  seems  to  be  preferable:  (a)  specimen  in 
asymmetric  four-point  bending  (Fig.6.11);  (b)  compact  specimen  for  three- 
point  loading  (Fig. 6. 15,  Case  3);  (c)  compact  specimen  for  two-point  loading 
(Fig. 6. 17).  These  specimens  provide  pure  Mode  II  at  the  crack  tip.  Some 
experimental  data  on  determination  of  fracture  toughness  (Table  2.1) 
demonstrate  that  the  -values  for  metallic  materials  are  normally  higher 
than/Cjc  .  However,  the  results  of  Tohgo  et  al.  (1990)  show  that  it  is  not  a 
common  rule.  For  brittle  materials  (e.g.,  perspex,  AI2O3  ceramics),  the  ratio 
^nc/^ic  is  slightly  less  than  unity.  Note  also  that  no  unique  fracture 
criterion  can  be  employed  to  predict  the  crack  initiation  angle  for  brittle 
and  ductile  materials.  The  maximum  tensile  stress  or  minimum  strain 
energy  density  criteria  can  be  used  in  predicting  a  crack  path  in  brittle 
materials,  while  the  maximum  shear  stress  criterion  seems  to  control  crack 
initiation  in  metals. 
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Figure  Captions 


Figure  1.  Coordinate  systems  at  a  crack  tip. 

Figure  2.  Schematic  illustration  of  Irwin’s  correction  for  plasticity. 

Figure  3.  Plastic  zone  size  ahead  of  the  crack  tip  in  Mode  II  based  on  the 
Dugdale  model  and  Irwin’s  correction  for  plasticity. 

Figure  4.  Plastic  zone  shapes  at  the  crack  tip  in  Mode  II  based  on 

the  Huber-Mises  and  the  Tresca  yield  criteria. 

Figure  5.  Plastic  zone  shapes  at  the  crack  tip  in  Mode  II  based  on  the 

two-parameter  representation  of  the  stress  field  and  the 
Huber-Mises  yield  criterion. 

Figure  6.  Crack  geometries  in  Mode  II. 

Figure  7.  Calibration  curves  of  the  stress  intensity  factor  for  the  compact 
shear  specimen.  Fig. 6. 15  ia=Pt{tH),  t  is  the  specimen  thickness), 
(a)  K II  for  the  Case  2  of  loading  at  different  sizes  of  the 
specimen,  after  Hoyniak  and  Conway  (1979).  (b)  Kn  and  S  for 
different  loading  schemes,  after  Savruk  (1988). 

Figure  8.  HRR  stress  and  strain  fields  for  Mode  II  crack  in  plane  strain. 

Figure  9.  HRR  stress  and  strain  fields  for  Mode  II  crack  in  plane  stress. 

Figure  10,  Variation  of  the  I(N)'factor  versus  N  for  plane  strain  and  plane 
stress. 

Figure  11,  Stress-strain  curves  for  technically  pure  iron  in  shear  tests  at 
different  strain  rates  (data  of  Klepaczko,1969). 

Figure  12.  (a)  Rate  sensitivity  of  the  lower  yield  stress  in  shear  for  mild 
steel  and  technically  pure  iron; 

(b)  Rate  sensitivity  of  constants  in  the  Ramberg-Osgood 
equation  for  technically  pure  iron. 

Figure  13.  Stress-strain  curves  for  titanium  Ti6A14V  alloy  from 
compression  tests  at  different  strain  rates  (data  of  Maiden  and 
Green.  1966). 

Figure  14.  Rate  sensitivity  of  the  yield  stress  in  shear  and  material 

constants  in  the  Ramberg-Osgood  equation  for  Ti6A14V  alloy. 

Figure  15.  Stress-strain  curves  for  A1  6061-T6  alloy  from  compression 
tests  at  different  strain  rates  (after  Maiden  and  Green,  1966). 
Conversion  factor  for  stress:  1  ksi  =  6.895  MPa. 

Figure  16.  Rate  sensitivity  of  the  yield  stress  in  shear  for  Ai  6061-T6 

alloy  (data  of  Jiang  and  Chen,  1974). 

Figure  17.  Variation  of  parameters  which  characterize  the  geometry  of 

plastic  zone  near  the  crack  tip  versus  logarithm  of  strain  rate 
for  a  semi-infinite  crack  in  mild  steel:  (a)  Xappi/xys^O--'*  (b) 
^appl^"*’ Y  s  •  5  • 

Figure  18.  Variation  of  parameters  which  characterize  the  geometiy'  of 

plastic  zone  near  the  crack  tip  versus  logarithm  of  strain  rate 


for  a  finite  length  edge  crack  in  mild  steel:  (a)  Tappi/'^YS=0.2;  (b) 

^appl/^YS=0.5. 

Figure  19.  Variation  of  parameters  which  characterize  the  geometry  of 

plastic  zone  near  the  crack  tip  versus  logarithm  of  strain  rate 
for  a  semi-infinite  crack  in  Ti6A14V  alloy:  (a)  tappi/xYS=0-2;  (b) 

"^app/”^  YS~®  •  ^  • 

Figure  20.  Variation  of  parameters  which  characterize  the  geometry  of 

plastic  zone  near  the  crack  tip  versus  logarithm  of  strain  rate 
for  a  finite  length  edge  crack  in  Ti6A14V  alloy:  (a)  tappi/xYs=0-2; 
(b)  ^appl/'CYS=0-5- 

Figure  21.  Boundaries  of  equal  equivalent  plastic  strain  at  the  crack  tip  in 
Mode  II  in  mild  steel  versus  strain  rate:  (a)  ej*  =  0.03; 

(b)  £P  =  0.06. 

Figure  22.  Boundaries  of  equal  equivalent  plastic  strain  at  the  crack  tip  in 
Mode  II  in  Ti6A14V  alloy  versus  strain  rate,  eP  =  0.02. 

Figure  23.  Finite-element  solutions  on  variation  of  the  plastic  zone 

boundary  versus  strain  hardening  exponent  at  the  crack  tip  in 
Mode  II:  (a)  after  Shih  (1974);  (b)  after  Pan  and  Shih  (1990). 

Figure  24.  Finite-element  solution  on  variation  of  the  plastic  zone 

boundary  versus  load  level  at  the  crack  tip  in  Mode  II  (after 
Banks-Sills  and  Sherman,  1990):  (a)  x^ppi/xys  =  0.35;  (b) 
^appl^^YS  =  0.45;  (c)  Xappi/Xys  =  0-48;  (d)  Xappi/XYS  =  0.52. 

Figure  25.  Boundaries  of  equal  equivalent  stress  at  the  crack  tip  in  Mode 
n  for  0=3/7,  N=3  (a)  and  N=13  (b). 

Figure  26.  Boundaries  of  equal  equivalent  plastic  strain  at  the  crack  tip  in 
Mode  II  for  a=3/7,  N=3  (a)  and  N:=13  (b). 


Figure  1.  Coordinate  systems  at  a  crack  tip. 


Figure  2.  Schematic  illustration  of  Irwin  s  correction 
for  plasticity. 


—  Dugdale’s  model,  edge  crack  of  finite  length 
X  Dugdale's  model,  semi-infinite  crack 

—  Irwin's  correction  tor  plasticity 
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Figure  3.  Plastic  zone  size  ahead  of  the  crack  lip  in  Mode  11  based 
on  the  Dugdale  model  and  Irwin's  correction  for 
plasticity. 
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Figure  6.  Crack  geometries  in  Mode  11. 
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Figure  7.  Calibration  curves  of  the  stress  intensity  factor  for  the 
compact  shear  specimen.  Fig. 6. 15  {cs-P/itH),  t  is  the 
specimen  thickness),  (a)  ATjj  for  the  Case  2  of  loading  at 
different  sizes  of  the  specimen,  after  Hoyniak  and 
Conway  (1979).  (b)  Kn  and  S  for  different  loading 
schemes,  after  Savruk  (1988). 
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Figure  8.  HRR  stress  and  strain  fields  for  Mode  II  crack 
in  plane  strain 


Plane  stress  fields,  N  =  3. 
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Plane  stress  fields,  N  =  3 


Plane  stress  fields,  N  ^  5. 
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Plane  stress  fields,  N  =  5. 


Plane  stress  fields,  N  =  7. 

- 1 - 1 - \ - 1 - 

60  90  120  150 

e,  c 


E 

E 

E‘ 


Plane  stress  fields,  N  =  7. 
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Plane  stress  fields,  N  =  10. 
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Figure  10.  Variation  of  the  I(N)’factor  versus  N  for  plane  strain  and 
I  plane  stress. 
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Figure  11.  Stress-strain  curves  for  technically  pure  iron  in  shear 
tests  at  different  strain  rates  (data  of  Klepaczko,1969). 


IS 


O  Campbell  sno  Ferguson  il97C),  0.12%  C 
©  T anaka  and  Kincsmta  ( 1 S67),  0.06%  C 
©  Klepaczko  (1963),  tecnnically  pure  iron 


(a)  log-/  [s' I 


(b)  log  Y  [s  1 

Figure  12.  (a)  Rate  sensitivity  of  the  lower  yield  stress  in  shear  for 
mild  steel  and  technically  pure  iron; 

(b)  Rate  sensitivity  of  constants  in  the  Ramberg- 
Osgood  equation  for  technically  pure  iron. 
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Stress-strain  curves  for  titanium  Ti6A14V  alloy 
from  compression  tests  at  different  strain  rates  (data  of 
Maiden  and  Green,  1966). 
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Figure  14.  Rate  sensitivity  of  the  yield  stress  in  shear  and 
material  constants  in  the  Ramberg-Osgood  equation  for 
Ti6A14V  alloy. 
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Figure  15.  Stress-strain  curves  for  A1  6061-T6  alloy  from 
compression  tests  at  different  strain  rates  (after  Maiden 
and  Green,  1966).  Conversion  factor  for  stress:  1  ksi  = 
6.895  MPa. 
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Figure  17.  Variation  of  parameters  which  characterize  the  geometry 
of  plastic  zone  near  the  crack  tip  versus  logarithm  of 
strain  rate  for  a  semi-infinite  crack  in  mild  steel:  (a) 
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Figure  18a 
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Figure  18.  Variation  of  parameters  which  characterize  the  geometry 
of  plastic  zone  near  the  crack  tip  versus  logarithm  of 
strain  rate  for  a  finite  length  edge  crack  in  mild  steel;  (a) 
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Figure  19.  Variation  of  parameters  which  characterize  the  geometry 
of  plastic  zone  near  the  crack  tip  versus  logarithm  of 
strain  rate  for  a  semi-infinite  crack  in  Ti6Ai4V  alloy;  (a) 
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Figure  20b 


Figure  20.  Variation  of  parameters  which  characterize  the  geometry 
of  plastic  zone  near  the  crack  tip  versus  logarithm  of 
strain  rate  for  a  finite  length  edge  crack  in  Ti6A14V  alloy; 

tappi/tYS^O-^'i  (b)  'tappi/tYS— 


Figure  22.  Boundaries  of  equal  equivalent  plastic  strain  at  the  crack 
tip  in  Mode  II  in  Ti6AI4V  alloy  versus  strain  rate,  ej  = 
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Figure  23.  Finite-element  solutions  on  variation  of  the  plastic 
zone  boundary  versus  strain  hardening  exponent  at  the 
crack  tip  in  Mode  II:  (a)  after  Shih  (1974);  (b)  after  Pan 
and  Shih  (1990). 
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Figure  24.  Finite-element  solution  on  variation  of  the  plastic 
zone  boundary  versus  load  level  at  the  crack  tip  in  Mode 
II  (after  Banks-Sills  and  Sherman,  1990);  (a)  tappi^^vs  = 
0.35;  (b)  XappAvs  =  0-^5;  (c)  tappi/xys  =  0.48;  (d)  Xappi/xys  = 
0.52. 


Figure  26.  Boundaries  of  equal  equivalent  plastic  strain  at  the  crack 
tip  in  Mode  II  for  a=3/7,  N=3  (a)  and  N=13  (b). 


